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Abstract. We present an updated review of Lorentz invariance tests in Effective 
field theories (EFT) in the matter as well as in the gravity sector. After a general 

CO . discussion of the role of Lorentz invariance and a derivation of its transformations along 

the so called von Ignatovski theorem, we present the dynamical frameworks developed 
within local EFT and the available constraints on the parameters governing the Lorentz 
breaking effects. In the end, we discuss two specific examples, the OPERA "affaire" 
and the case of Hof ava-Lifshitz gravity. The first case will serve as an example, and a 
f^j\. caveat, of the practical application of the general techniques developed for constraining 

Lorentz invariance violation (LIV) to a direct observation potentially showing these 

effects. The second case will show how the application of the same techniques to a 

^ , specific quantum gravity scenario has far fetching implications not foreseeable in a 

•^ ' purely phenomenological EFT approach. 
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1. Introduction 

Quantum gravity lias been a frustrating endeavour for more than sixty years and in 
spite of many technical progresses it is a program far from completion. We have 
nowadays many approaches to quantum gravity, string theory, loop quantum gravity, 
causal dynamical triangulations, causal sets, quantum graphity, just to name a few, 
but these are radically different approaches, driven by underlying different guiding 
ideas and even general lessons about the nature of spacetime beyond the Planck scale 
are difficult to draw. Why so? The answers probably lies in the fact that Quantum 
Gravity has been a collective effort in theoretical physics like no other before, a search 
without lighthouses i.e. without observations to guide us. Indeed, we expect QG 
effects at experimentally /observationally accessible energies to be extremely small, due 
to suppression by the Planck scale Mpi = ^Jhc/G^ ^ 1.22 x 10^^ GeV/c^. In this sense 
it has been considered (and it is still considered by many) that only ultra-high-precision 
(or Planck scale energy) experiments would be able to test quantum gravity models. 

Nonetheless, no serious theoretical physicist or cosmologist can nowadays be 
satisfied of our current theoretical toolbox. Indeed, the standard model of particle 
interactions (SM) and classical General Relativity (GR) are extremely powerful tools 
and are rightly listed among the most impressing achievements of mankind. Yet, these 
theories are not enough for answering many pressing questions. Our observations are 
starting to probe physics just 10'^^ seconds after the Big Bang asking for theories able 
to set the initial conditions of our Universe. Furthermore, our cosmology makes sense 
within our current theoretical framework only at the price of a very preposterous choice 
of matter-energy content, we do understand only about 4% it. Even worse GR fails to be 
a predictive theory in many interesting and relevant situations. Indeed, many solutions 
of Einstein's equations are singular in some region, not least the very beginning of our 
Universe or the interior of black holes (for which we do have nowadays a solid evidence 
of existence). Moreover, there are honest classical solutions of the Einstein equations 
that contain closed time-like curves, which would allow traveling back and forth in time 
with the associated causal paradoxes. Finally, the problem of black-hole evaporation 
considered just within the framework of semi-classical gravity clashes with quantum 
mechanical unitary evolution. 

The conclusion is one and the same: we need nowadays much more than in the past 
a theory of quantum gravity, something able to describe the nature and the dynamics 
of spacetime beyond the Planck scale. So, even more than in the past any attempt 
to put boundaries and constraints on our intuitions/ideas/speculations about what lies 
behind the Planck scale are necessary and should not be dismissed as simple minded or 
premature. 

The broad and multi-faceted field of "quantum gravity phenomenology" has 
represented in the last twenty years a difficult, controversial and painstaking attempt to 
do this. To grasp whatever could be grasped by observations and to tests any (more or 
less crazy) idea could be advanced and tested about Planck scale effects and their low 
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energy imprint on reality. 

In fact, quantum gravitational models beyond GR have shown that there can be 
several of what we term low energy "relic signatures" of quantum gravitational effects 
which would lead to deviation from the standard theory predictions (SM plus GR) in 
specific regimes. Some of these new phenomena, which comprise what is often termed 
"QG phenomenology" , include tests of 

• Quantum decoherence and state collapse [1] 

• QG imprint on initial cosmological perturbations [2] 

• Cosmological variation of couplings [3l H] 

• TeV Black Holes, related to extra-dimensions [5] 

• Planck scale spacetime fuzziness [6] 

• Generalized uncertainty principle [3 [HI E] 

• Violation of discrete symmetries [10] 

• Violation of space-time symmetries [11] 

This review will focus on a special subset of the last item i.e. Lorentz violations. 

1.1. Why testing Lorentz invariance? 

Lorentz invariance lies at the roots of both the SM and GR and as we shall see soon it is 
really an outcome of our deepest understanding of the nature of space and time. So why 
we should sacrifice it on the altar of quantum gravity, isn't this a too high price to pay? 
The answer to this question is twofold. First, there is a practical issue: exactly because 
Lorentz Invariance (LI) is a cornerstone of our current understanding of reality we do 
have a compelling duty to test it as far as we can. Secondly there are good reasons 
to investigate about the compatibility of this cherished symmetry with Planck scale 
physics. In fact, many models of quantum gravity involve some form of discretisation 
of spacetime, something that a priori is hard to reconcile with e.g. boost invariance. Of 
course, one can envisage a discretisation of spacetime which is by construction LI, like 
e.g. in causal sets. However Lorentz invariance of the sub-Planckian structure of a QG 
model does not per se guarantee the preservation Lorentz invariance at intermediate 
scales. Let us dwell a bit more carefully on this point. 

All our models of QG have to solve a crucial issue related to the emergence of 
a classical spacetime at sufficiently low scales. It is often assumed that as far as the 
underlying structure does not involved a preferred frame so the emergent structures 
will. This is however not supported by simple toy models that one can construct in this 
sense. To cite just a simple example, it was shown within the so called Analogue Gravity 
framework [12l [13] that a system of relativistic bosons that undergoes a Bose-Einstein 
condensation is characterised by low energy collective excitations (quasi-particles) that 
experience the condensate as an effective spacetime and are characterised by a relativistic 
dispersion relation with the speed of light replaced by the speed of sound in the BEG, 
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w^ = c^/c^. However, the same excitations at high energy have cumbersome dispersion 
relations that break what we might call this "acoustic" LI and interpolate finally with 
a UV fully relativistic dispersion relation this time characterised by the speed of light 
as the limit speed. While this toy model cannot be accounted in any way as a model of 
quantum gravity (as it shows an emergent spacetime but cannot reproduce and emergent 
GR dynamics |j) it is definitely making a point about the fact that an underlying Lorentz 
invariance of the system is not sufficient (and could be not even necessary) in order to 
provide an almost exact Lorentz invariant physics in the emergent spacetime. 

2. The von Ingnatovski theorem 

In discussing the viable test of Lorentz invariance it might be fruitful to take some time to 
express first what we are actually probing when we cast the aforementioned constraints. 
It is in fact often associated with LI the observation that the speed of massless particles 
(i.e. the speed of light) is the same independently of the chosen reference system in 
which the measurement is performed. 

This is true and false at the same time. While it is surely true that a route to the 
derivation of Lorentz transformations is to assume a priori the invariance of the speed 
of light, it was also soon understood [171 CHI CHI 1201 E] that no such assumption was 
needed and that indeed the basic transformations of Special Relativity could be derived 
starting from basic assumptions about the fundamental nature of spacetime. 

In this sense a first step consists in realising that no metaphysical, abstract, meaning 
needs to be associated to spacetime. Indeed, this can be defined operationally by 
replacing the abstract notions of "time" and "space" with the physical notions of 
"duration" (of a physical phenomenon as measured by a clock) and distance (e.g. of 
a physical object as measured by a suitable rod). It should be clear that this 
"concretisation" of spacetime implies that its observed features and symmetries will be 
always related to properties and symmetries of the clocks and rods used to determine 
durations and distances and of the measurement protocols adopted. The introduction 
of a set of observers endowed with standard clocks and rods and the agreement on a 
suitable procedure of synchronisation of such clocks allows to rigorously define such a 
spacetime and observations tells us that locally (when gravity can be safely neglected) 
space is isotropic and space-time are homogeneous. Typically the notion of time as 
measured by our clocks is taken for simplicity in such a way that the motion of an 
isolated particle is uniform x = 0. A system of observers satisfying this request is said 
to be inertial. One can then prove that if /C is an inertial system of reference then any 
other system K, moving with respect to it of uniform rectilinear motion is inertial. 

The above premises lay down the basic elements for an axiomatic derivation of the 
Lorentz transformations [22l |23l [2l] , here we shall discuss it following the logic proposed 
in [25]. Let's summarise the necessary assumptions. 



I Sec however |14| 1151 I16j for anafogue gravity or an anafogue gravity inspired models showing possible 
emergence of gravitational-like dynamics in these systems. 
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• Spatial and temporal homogeneity 

• Spatial isotropy 

• The Relativity Principle (equivalence of inertial frames) 

• Pre-Causality (the time-ordering of two events along an observer worldline does not 
change in different systems of reference) 

Let us start by considering two inertial frames ]C{0,x,y,z,t) and }C{0,x,y,z,t) 
where O and O are respectively the origins of the two frames and we shall assume that 
/C moves w.r.t. K, with a uniform rectilinear speed v (which, for simplicity, we can take 
along the x direction). Note that a priori one can only presume that v G (— c_, c+) = J 
with c_,c_|_ > but possibly infinite. Similarly one can say that fC moves in /C with 
velocity v (for the moment let us not assume anything about the v,v relation which we 
shall discuss later). 

We can always choose the two origins to coincide at t = t = i.e. t = x = y = z = 
iii=x = y = z = and we shall derive now what are the transformation laws relating 
the X and x coordinates, x^ = f^(t, x, y, z, v). 

Now homogeneity in space and time requires that the above mentioned 
transformation laws are linear. This can be easily understood via the following 
reasoning: homogeneity implies that time or space intervals measured say in K, can 
depend only on the corresponding intervals in /C, not on the precise instants and 
positions when and where the measurements are made. Hence 

dx^ = — — dx'^ with — — = const . (1) 

dx" dx" ^ ' 

I.e. /^ must be a linear function of the coordinates x^ . Specialising for simplicity to 

the 1+1 case, the above statement translates in the simple form for the (t, x) — )■ (t, x) 

transformations 

t = A{v)t + B{v)x , 

x = C{v)t + D{v)x . (2) 

Let us now ask how the motion of the origin O appears in the system of reference 
]C. We chose coordinates so that at t = the origin O had to be in x = 0, hence at any 
later time xq = vi. However, we do also know that at any time t, O has to be located 
at X = hence from ([2]) we get for the coordinates of an event that takes place at O 

i = A{v)t , 

xo = C{v)t , (3) 

from which (using xq = vt) easily follows C{v) = vA{v). 

We can now repeat the same argument in describing the motion of O w.r.t. /C. This 
will move with speed v in K, and reasoning as before we have Xq = vt and given that O 
is always located in x = in ^ from ([2]) we get 

C{v)t + D{v)xo = 0, (4) 
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from which (using Xq = vt) easily follows C{v) = —vD{v) which, given the previous 
result, implies 

D[v) = - ^ = -Ia{v) . (5) 

We can then express the general transformation relating the two systems of reference 
as 





(6) 
where A(f ) is the matrix 

A(«) := Mv) ( ' «1 ] , (7) 



and we have introduced the ratio C,{v) '■= B{v)/A{v). 

We can now make use of the isotropy of space (see [25] for a more general derivation 
not assuming a priori isotropy) to relate v and the reciprocal velocity v. These are 
obviously related by the Relativity Principle which implies that ii v = F{v) then 
V = F{v) or equivalently F = F~^. But this fact per se does not imply v = —v, 
the so called Reciprocity Principle, as is normally postulated. (Note that actually any 
F symmetric under diagonal reflections x — ?■ —x would be Ok.) Indeed, one needs in 
addition to make explicit use of the isotropy of space requiring the symmetry of the 
transformation under an inversion of the spatial axis x [26l [27] . (It is interesting to note 
that von Ignatovski erroneously thought that the reciprocity principle can be derived as 
a consequence of the Relativity Principle.) 

Hence, in conclusion, Eq. ([7]) further simplify into 



A{v):=Aiv) ; ^V^^ , (8) 




Obviously, A(0) must be the identity matrix, so A{0) = 1 and ^(0) = 0. We can 
also check that invariance of the transformation under spatial inversion (e.g. x — ?• —x, 
X — )■ —X, u — 7- —u) or time reversal requires A{v) = A{—v) and ^{v) = —^{—v). 

We are now ready to enforce the group structure into our set of transformation by 
appealing to the Relativity Principle. The way to do so is to introduce a third system 
of reference K. which moves with constant velocity u w.r.t. /C and u w.r.t. K.. We expect 
u to be given by some composition law 

u = $(M,f) , (9) 

while the transformation relating K. to K, will be 

(10) 
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This, together with Eq. ([B]), also tells us that the transformation from K, to K, will have 
the form 

'*j=AW-'A(.)(*j. (11) 

However, we also know that the transformation 

:)-A(«)(j). (12) 

must independently hold which is tantamount to say that 

A{v)'^A{u) = A{u) =^ A{v)A{u) = A{<l>{u,v)) . (13) 

This equation implies the group structure for our set of coordinate transformations. 
Note that if u* is some value of the velocity u so that u = ^{u*,v) = then 
A{v)A{u*) =A(0) = I 

Let us now differentiate w.r.t. v equation (IT3|l and calculate it in f = 0. We get 



A'(0)A(u) = ip(u)A'(u) where ip(u) 



dv 
It is easy to see that this equation admits the solution 



(14) 



v=0 



A(n) = e'^(")^'(°) where h{u) = [ -^ . (15) 

Jo ^W 

Now we can use Eq. ([8]) to compute 

A'(0) = A'{0)I + M where M = j ° ^'^°^ j . (16) 

The matrix M has the following algebraic properties: For any natural number k, 

Wf'' = (e'(0))^ I , Wf^^^ = (e'(0))^^ M . (17) 

Replacing these into Eq. (II 5p we find 



where: 



A(m) = e^(^)^'(°) ( ''^""^ ^'(0)^2(n) 

S2{u) Si{u) 



+ 00 ^ 



fc=0 



{2k)\ 



+ 00 , 

fc=0 ^ ~^ )■ 
We can now compare Eqs. ([8]) (with v replaced by u) and (TT8l) to find: 

A(M)=e'^(")^'(°)si(M); (21) 

«-^; (22) 

Si{u) 
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^iu) = e{0)u; (23) 

Note that Eq. ( I22l) contains implicitly the link between u and h{u), because its right- 
hand side depends on u only through h{u). Furthermore, we previously saw that 
A(u) = A{—u) due to spatial or time reversal invariance. However, the function h{u) is 
odd under this operation, h{—u) = —h{u), while Si{u) is evidently even. By consistency 
one then must ask A'{0) = which in the end implies 

Replacing Eq. ( l24l) and the corresponding expressions for v and $(m,w) into Eq. 
13]), we find: 



h{^{u,v)) = h{u) + h{v); (25) 

Si($(n,i;)) = Si(n)si(tO(l + r(0)) ; (26) 

and most importantly 

Equation fl27|) implies (p{u) = —^'{0) n^ + 1, furthermore tells us that that the condition 
^'(0) > must be satisfied in order to guarantee that $(m,w) does not diverge for any 
finite value of m or w. 

Finally, using the relationship s'2{u) = Si{u)/ip{u) (which follows from Eqs. ( ITSll . 
(dSD, ([20]) and ^) together with Eqs. ([22]) and <^{u) = -^'{0) u^ + 1, we obtain 

s'2(u) /1\ 1 1 ^'{u) 



S2(n) \u) -i'(Si)v?^-\ u 2(^(n) ■ ^ ' 

This can be immediately integrated to obtain, in a neighbourhood of the origin, 
■52 (it) = ±m/a/v5(m); hence S\{u) = ±/y^^p(u) and 

^^""^ ^ T^M ^ Vl - e'(0)n2 ' (29) 

where a multiplicative constant and the sign ambiguity have been fixed by the condition 
A{0) = 1. 

All we need now is to impose what we called the pre-causality axiom which formally 
translates into the condition di/dt > 0. In our set of transformations ([6]) this implies 
A{u) > \/u & J. However, we already know that ^'(0) > 0, hence for A{u) to be 
positive and real one has to require ^'(O)m^ < 1. 

We can then define a speed scale C^ = l/.^'(0) and require C^ > u"^ for any u and 
we can now write 

7(«) = A{u) = \ . (30) 

The speed scale C has two crucial properties 
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(i) From (127|) we can easily see that it is a fixed point in the composition law 
$(C,t,) = $(u,C) =C. 

(ii) The transformations ([6]) holds only for \u\ < C so that J = (— C, C). 

Note that we have not proved at all that signals cannot move faster than C\ In principle 
these particles/signal are not incompatible with Lorentz transformations and they do 
not a priori lead to causality violations |28] . 

Finally our transformations ([6]) will then take the familiar form 

A{v) := ^{v) ( I ''^f j , (31) 

It is only empirical observation which will tell us in the end if C is finite (Lorentzian 
relativity) or infinite (Galilean relativity). Furthermore, nowhere in our derivation of 
the Lorentz transformations we needed to appeal to the universality of the speed of 
light. Even worse we do not know a priori that the limit speed C has to coincide with 
such a physical entity. Let us stress that the determination of C is strictly related 
to the physics underlying our system of clocks and rods. This was for example nicely 
described in a pedagogical novel by Trautman [29] where bats using only sound wave 
to describe their spacetime would easily end up deriving the speed of sound as such an 
invariant /limit speed (see also [30] for an acoustic gravity analog concrete realisation of 
the same idea). 

In this sense the very notion of Lorentz invariance is detached by an abstract 
geometric entity as the Minkowski spacetime and firmly attached to a symmetry of 
the physics underling our tools for operationally define the spacetime. Conversely 
the Minkowski spacetime can be seen as a mathematical tool that summarises such 
properties. In this sense tests of departures from Lorentz invariance can of course be 
seen on the one side as probes of the microscopic structure of space-time but can as well 
be considered high precision tests of the physics underlying our most precise clocks and 
rods. 

3. A historical overview 

The previous discussion should have convinced the reader about the fundamental nature 
and insight provided by Lorentz invariance in modern physics. It is then not surprising 
that questioning and investigations regarding such fundamental symmetry appeared 
relatively early in last century literature. It is however only in recent times that 
speculations could be translated in full-fledged phenomenological studies. We shall 
here present an incomplete review of these developments. 

The possibility that Lorentz invariance violation (LIV) could play a role again in 
physics dates back by at least sixty years [311 |32l [331 Ell ESI ES] and in the seventies 
and eighties there was already a well established literature investigating how LI could 
be established at low energies without being an exact symmetry at all scales (see 
e.g. [3?! EHl EHl Hni HH 112] )• This stream of paper had a lasting influence but somehow 
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remained as a set of isolated attempts to go beyond the standard framework. In general, 
the gravitation theory community did not pay very much attention to the subject given 
the general expectation at the time that possible breakdowns of standard physics due 
to quantum gravitational effects were only to appear in particle interactions at energies 
of the order Planck mass Mpi. It was only in the nineties that it was clearly realised 
that there are special situations in which new effects could manifest also at lower energy. 
These situations were termed "Windows on Quantum Gravity". 

While the Planck mass stands off from any standard model or observational scale 
(Ultra High Energy Cosmic Rays, see, e.g., |l3l IB]) have E < lO^^ GeV ~ IQ-^Mp;, 
sometimes even tiny deviations from standard physics can be magnified into a significant 
effect when dealing with threshold reactions at high energies (but still well below the 
Planck scale), long distances of signal propagation, or other very sensitive phenomena 
(see, e.g., [H] for an extensive review). A partial list of these windows on (^G^ includes: 

• sidereal variation of LIV couplings as the lab moves with respect to a preferred 
frame or direction 

• cumulative effects: long baseline dispersion and vacuum birefringence (e.g. of signals 
from gamma ray bursts, active galactic nuclei, pulsars) 

• anomalous (normally forbidden) threshold reactions allowed by LIV terms 
(e.g. photon decay, vacuum Cherenkov effect) 

• shifting of existing threshold reactions (e.g. photon annihilation from Blazars, ultra 
high energy protons pion production) 

• LIV induced decays not characterised by a threshold (e.g. decay of a particle from 
one helicity to the other or photon splitting) 

• maximum velocity (e.g. synchrotron peak from supernova remnants) 

• dynamical effects of LIV background fields (e.g. gravitational coupling and 
additional wave modes) 

As usual in science, the realisation of this potentiality was a community effort, 
however one can identify some seminal papers that contributed to the take off of this 
field of research. Among these papers one can cite |15] that already in 1989 envisaged, 
within a string field theory framework, the possibility of non-zero vacuum expectation 
values (VEV) for some Lorentz breaking operators. This work led later on to the 
development of a systematic extension of the SM (that was later on called "minimal 
standard model extension" (mSME)) incorporating all possible Lorentz breaking, gauge 
invariant, power counting renormalizable (i.e. of mass dimension < 4) [36]. This 
provided a framework for computing in effective field theory the observable consequences 
for many experiments and led to much experimental work setting limits on the LIV 
parameters in the Lagrangian (see e.g. [47]). 

Meanwhile, other theoretical investigations provided new motivations for Lorentz 
breaking searches and constraints. Indeed, specific hints of LIV arose from various 
approaches to Quantum Gravity. Among the many examples are the above mentioned 
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string theory tensor VEVs ^5\ and space-time foam models |1S1 1121 EOl ED IS2] • It was 
this stream of research that finally lead to the realisation that constraints could be cast 
also on high energy violations of Lorentz invariance in the photon dispersion relation, 
using the aforementioned propagation over cosmological distance of light from remote 
astrophysical sources like gamma ray bursters (GRBs) and active galactic nuclei (AGN), 
was instead highlighted in |19] The field of phenomenological constraints on quantum 
gravity induced LIV was born. 

In fact, starting from this early work a stream of theoretical and phenomenological 
papers followed. For example, the influential papers by Coleman and Glashow 
[531 IHU [55] which brought the subject of systematic tests of Lorentz violation to the 
attention of of the broader community of particle physicists. Let me stress, that several 
other papers appeared in the same period, some of them anticipating many important 
results, see e.g. [561 EZ]- However at the time of their appearance were hardly noticed 
(probably seen by many as too "exotic"). 

In the very same period an important part of the revived interest LIV was due 
to the development of Analogue Gravity [58]. While the field was initiated in 1981 by 
Bill Unruh [59] it was only in the early 2000 that it was considered more systematically 
as an inspirational source for possible new effects associated to an "emergent gravity" 
framework. In particular, modified dispersion relations of the sort suggested in these 
analogue systems provided an ideal setting for possible testable signature on the cosmic 
microwave background of Planck physics [60l [6T1 [62] . 

The above field was joined by many other suggestions from novel QG 
models/calculations. For example, semiclassical spin-network calculations in Loop 
QG [63], non-commutative geometry [Ml [65l [66], some brane- world backgrounds [67] . 
Furthermore, during the last three decades there were several attempts to formulate 
alternative theories of gravitation incorporating some form of Lorentz breaking, from 
early studies [HHl [691 EOl UH [12] to full-fledged theories such as the Einstein- Aether 
theory [73l [ZH [75] and Hofava-Lifshitz gravity [761 EH EH] (which in some limit can be 
seen as a UV completion of the Einstein-Aether framework [79]). 

Coming back to the phenomenology side of this story, it was again in the early 
years 2000 that several works made an effort to provide a more systematic study both 
of LIV frameworks as well as of the available constraints (see e.g. [801 [8T1 [82]). In this 
sense a crucial contribution was the development of an effective field theory approach 
also for higher order (mass dimension greater than four), naively non-power counting 
renormalizable, operators. [§| This was firstly done for dimension 5 operators in QED 
[85] to be later on extended to dimension 6 operators [86] . 

Noticeably, in the very same period was advanced the idea that new relativity 
groups might exists while being able to incorporate an invariant length/energy scale 



§ Anisotropic scaling |83[ 1761 I84j techniques were recently recognised to be the most appropriate way of 
handling higher order operators in Lorentz breaking theories and in this case the highest order operators 
arc indeed crucial in making the theory power counting renormalizable. This is why we shall adopt 
sometime the expression "naively non renormalizable" . 
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to be identified with the Planck one. This is the stream of research associated to the 
so called Doubly/Deformed Special Relativity [871 [Ml ESI |90]). This field lead, quite 
recently, to more bold proposals such as Relative Locality jHH [921 [931 EH [95] . We shall 
here preeminently discuss tests of Lorentz violations within an EFT framework but we 
shall come back on these alternative framework in the end of this work. 

It should be stressed that this burst of activity was not only spurred by the simple 
realisation that some phenomenological constraints on Planck physics could be cast. 
It is in fact important to note that the very same period was characterised by a 
growing plethora of observational puzzles seemingly hinting towards new physics. For 
example, in cosmology these are the years of the striking realization that our universe 
is undergoing an accelerated expansion phase [961 ISZ] which apparently requires a new 
exotic cosmological fluid, called dark energy, which violates the strong energy condition 
(to be added to the already well known, and still mysterious, dark matter component). 

Also in the same period high energy astrophysics provided some new puzzles, first 
with the apparent absence of the Greisen-Zatsepin-Kuz'min (GZK) cut off [981 [HH] (a- 
suppression of the high-energy tail of the UHECR spectrum due to UHECR interaction 
with CMB photons, see section 17.31 below for further discussion) as claimed by the 
Japanese experiment AGASA [100] . later on with the so called TeV-gamma rays crisis, 
i.e. the apparent detection of a reduced absorption of TeV gamma rays emitted by 
AGN |101] . Both these "crises" later on seemed to subside or at least alternative, more 
orthodox, explanations for them were advanced. However, their undoubtedly boosted 
the research in the field at that time and recent observational developments leave these 
two puzzles still open to further scrutiny. 

In this sense we cannot omit here the recent "Opera affaire", i.e. the false detection 
by the CERN-LNGS based experiment OPERA |102] of superluminal propagation of 
muonic neutrinos and Lorentz EFT (see e.g. ^U3\ [IDl MS\ ^Q6\ [Ml [I08]). While it 
is nowadays clear that the measurement was flawed due to unaccounted experimental 
errors, it did play an important role in propelling further activity in Lorentz breaking 
phenomenology which provided useful insights for future searches. We shall come back 
on this issue later using this case as a test field of the general techniques developed 
within a EFT framework. 

4. Matter Lorentz breaking effective field theory 

Let us now discuss in detail the specific frameworks introduced within the EFT approach 
to Lorentz breaking phenomenology, which in the case of the matter sector is sometime 
referred to as the Standard Model Extension (SME). This consists in supplementing the 
standard model of particle physics with all the possible Lorentz violating operators that 
can be written without changing the field content or violating gauge symmetry. The 
operators appearing in the SME can be conveniently classified according to their mass 
dimension and behaviour under CPT. Note that Lorentz violation does not imply CPT 
violation for local EFTs |109] . while CPT violation does imply Lorentz violation in local 
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EFT [TTU] . 

Within the SME it is often considered the specific subset of the rotationally 
invariant LIV operators. One reason for restricting our analysis to this specific subgroup 
of operators is of course simphcity and clarity. In the full rotation breaking theory there 
appears very many independent (and often dimensionfuU) coefficients which in turn 
lead to long tables of constraints from which it is sometime difficult to extract a sound 
physical intuition. Also theoretical arguments for LIV generally imply that it may arise 
in QG from the presence of a short distance cutoff which suggests just breaking of boost 
invariance, with a preferred rest frame||jl Furthermore, it would be very difficult for 
a theory which breaks rotation invariance to preserve boost invariance - if a spacelike 
four-vector is introduced to break rotation invariance, the four-vector also breaks boost 
invariance. Hence, a constraint on pure boost violation is, barring a conspiracy, also 
a constraint on boost plus rotation violation. We shall then from here on simplify our 
EFT treatment by always specialising to a rotational invariant setting. 

Once rotational invariance is assumed, one can easily see that all LIV tensors must 
reduce to suitable products of a time-like vector field, usually denoted m° which is usually 
taken to be unit (so that in the frame of the observer whose world line is tangent to m", 
the latter has components (1,0,0,0)). This allows us to express constraints solely in 
terms of the numerical coefficient involved in any M°-matter interaction term. Of course, 
the actual direction of u" is technically arbitrary. However, the common choice, which 
we make here, is to define u" to be aligned with the rest frame of the cosmic microwave 
background. Since the boost factor of an Earth centred frame is only ~ 10~^ |114] 
with respect to the CMB frame, any rotation breaking effects generated by the Earth's 
relative motion with respect to the CMB will be naturally suppressed by a factor of 
~ 10~^ relative to rotationally invariant effects and so we ignore them. 

4.1. mSME 

The Minimal Standard Model Extension or mSME, is the subset of the SME which 

considers only mass dimension 3 and 4 operators (naively power counting renormalizable 

operators). For fermions and photons the rotationally invariant LIV operators then take 

the form 

- - 1_^1_^ 

- au^ij'j^^ -bu^^p'j^Y^ + l^^cUf.u^ij'j'' D" i^+'-idu^u^^'j^'j^ D" ^j , (32) 

- ^{kF)u^vx>.'>^,F^^F^^'' + h,AFU^e^^p,A^FP\ (33) 

The dimension 3, CPT odd kAp term generates an instability in the theory and so 
is generally set to zero. The a term can be absorbed by shifting the phase of the fermion 
field and so we will ignore it temporarily (note however that it can be relevant in gravity, 
as we shall see later, and in the case of flavour dependent LIV in neutrinos |115j ). 

II See however |111) for an example where (coarse grained) boost invariance is preserved in a discrete 
model, and [112| I113J for a study of how discreteness may be compatible with Lorentz symmetry in a 
quantum setting. 
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The corresponding high energy (Mpi ^ E ^ m) dispersion relations for QED can 
be expressed as (see e.g. [11] and references therein) 

El = ml+p^ + jf^p + fi^^p^ electrons (34) 

^^ = (1 + 42)^ photons (35) 

where /i = — 2fo, /i = — (c — ds), and /^ = kp/2 with s = ±1 the helicity state 
of the fermion [11]. The anti-fermion dispersion relation is the same as (!34|) with the 
replacement p — )■ — p, which will change only the fe term. 

Note that the typical energy at which new phenomenology should start to appear 
is quite low. In fact, taking for example fe ~ 0(1), one finds that the corresponding 
extra-term is comparable to the particle mass m precisely at p ~ ?7i. Even worse, for 
the linear modification to the dispersion relation, we would have, in the case in which 
/i^^ ~ 0(1), that pth ~ mVMpi ~ lO'^'^ eV for electrons. 

4-2. Rotational invariant CPT odd dimension 5 ops 

An alternative approach within EFT is to study naively non-renormalizable operators. 
It is nowadays widely accepted that the SM could just be an effective field theory 
and in this sense its renormalizability is seen as a consequence of neglecting some 
higher order operators which are suppressed by some appropriate mass scale (generally 
the Higgs or SUSY one). One can then speculate that at least some of these non- 
renormalizable operators can be generated by quantum gravity effects (and hence be 
naturally suppressed by the Planck mass), being associated to the violation of some 
fundamental space-time symmetry like local Lorentz invariance. 

In [85] it was found that there are essentially only three operators of dimension five, 
quadratic in the fields, that can be added to the QED Lagrangian preserving rotation 
and gauge invariance, but breaking local LI. I5J 

These extra-terms, which result in a contribution of 0{E/Mpi) to the dispersion 
relation of the particles, are the following: 

^ -u'^FUu ■ d){ur,Fn + -i-n™^7„(Ci + C,^,){u ■ 9) V , (36) 



2Mpi ' '' ' 2Mpi 

where F is the dual of F and .^, Ci,2 are dimensionless parameters. All these terms also 
violate the CPT symmetry. More recently, this construction has been extended to the 
whole SM [TTB] . 

From fl36|l the dispersion relations of the fields are modified as follows. For the 
photon one has 

ul = e± -^e , (37) 

Mpi 

% Actually these criteria allow the addition of other (CPT even) terms, but these would not lead to 
modified dispersion relations (they can be thought of as extra, Planck suppressed, interaction terms) 
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(the + and — signs denote right and left circular polarisation), while for the fermion 
(with the + and — signs now denoting positive and negative helicity states) 

El=p' + m' + v^^, (38) 

with r]± = 2((^i±C2)- For the antifermion, it can be shown by simple "hole interpretation" 
arguments that the same dispersion relation holds, with rjg^ = —rjq^ where q and q 
denote respectively anti-fermion and fermion [821 1117] . 

As we shall see, observations involving very high energies can potentially cast 0(1) 
or stronger constraint on the coefficients defined above. A natural question arises then: 
what is the theoretically expected value of the LIV coefficients in the modified dispersion 
relations shown above? This problem could be further exacerbated by renormalization 
group effects, which could, in principle, strongly suppress the low-energy values of the 
LIV coefficients even if they are 0(1) at high energies. Let us, therefore, consider the 
evolution of the LIV parameters first. This issue was addressed in |116] by calculating 
the renormalization group equations for QED and the Standard Model extended with 
dimension-five, CPT odd, operators. The running was found to be only logarithmic and 
therefore low energy constraints are robust: 0(1) parameters at the Planck scale are 
still 0(1) at lower energy. Moreover, they also show that 77+ and ?7_ cannot, in general, 
be equal at all scales. Similar results were found also in the mSME and are expected at 
all orders. 

4.3. Rotational invariant CPT even dimension 5 and 6 ops 

The CPT even, rotationally invariant mass dimension five and six LIV terms has been 
computed in [86] using the same procedure described before for dimension 5 LIV. (This 
has been later extended to non-rotationally invariant operators in [118] 1115] .) 
With this ansatz one has for the photons 

- ^/3f i^"^^M«'^(^ ■ d)F^, . (39) 

While for fermions the mass dimension five and six operators are 

^ ii{u-Df{afPL + afPR)^ (40) 



Ml 



PI 



I 



^,-i^{u ■ D)\u . i){afPL + afPR)i^ (41) 



I —, 



Ml, 



A,(6)p , A,(6) 



i){u ■ D)U{u ■ iWI'Pl + «F^r)^ • (42) 



From these operators, the dispersion relations of fermions and photons can be easily 
computed, yielding 

E' - p' - m^ = "^j^iE + sp) + '^^{E - sp) + 
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c.^-fc^ = /3(«)-^, (44) 

where m is the electron mass and where s = cr • p/|p| is the hehcity of the fermions, 
and we have neglected terms of order m^/Mpj stemming from the a^^-* terms as they are 
highly suppressed. 

We can further simplify eq. ( 144|) by nothing that high-energy {E ~ p) fermions 
states are almost exactly chiral, and by grouping terms with same powers of the 
momentum. In this case we get 

E'=p' + m^ + ^rj^%^ + rj^:^^ (45) 

where R = +, L = — and we have labelled r]^'^^ is the dispersion coefficient of the LIV p^ 
term in the dispersion relation for the fermion. We choose r]^^^ as the fermion coefficient 
symbol as this nomenclature is common in the literature. Similarly, we shall use ^^"^ for 
the generic dispersion coefficient for a photon (so in (H4|) we shall take /J'-^^ = ^^^■'). 

Let us briefly comment on the so obtained dispersion relations. Eq. (H5ll 
distinguishes different fermion helicities albeit being generated via the addition of CPT 
even terms, this is due to the presence of LIV terms among those displayed in eq. fl39|) 
which are odd under P and T and hence can distinguish among helicity states. 

In (H5l) the quadratic modiflcation generated by the dimension flve operator is 
suppressed by a factor of order m/Mpi and hence it can be safely neglected, provided that 
E > y/mMp\. Note, however, that this term can be taken as an example of a dimension 
4 LIV term with a natural suppression, which for electron is of order rrie/Mpi ~ 10~^^. 
One should take this as a cautionary note as it shows clearly that for terms at this order, 
where no evident suppression is factorized out, the strength of a constraint should be 
judged only a posteriori, i.e. on the base of a model predicting the natural size of the 
dimensionless coefficient in front of the p^. In particular, to date the best constraint for 
a rotationally invariant electron LIV term of dimension 4 is O(10~^^) |119] which would 
not be enough to effectively constrain the above mentioned term, given that its natural 
size is 10"^^. 

By CPT, the dispersion relation of the anti-fermion is given by ( H^ . with the 
replacements s — )■ — s and p — ?■ —p. If q, q denote a charge fermion and anti-fermion, 
then the relevant anti-fermion coefficient r]^ is such that rn = r]q-^ , where q^ indicates 
an anti-fermion of positive/negative helicity (and similarly for the q±). This should be 
compared with what we got from a standard "hole interpretation" argument in the case 
of CPT odd dimension flve operators. It is easy to realise that this can be generalised 
via the same argument to conclude that for arbitrary "n" order terms [n being the 
power of the momentum) one would expect r]^ = (— l)"?7g" . This different behaviour 
between even and odd powers of "n" type dispersion relations leads to quite distinct 
phenomenologies as we shall see later. 
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4.4- The Neutrino sector 

The generic neutrino LIV operators, at any mass dimension, have been categorized in 
[115] . Also in this case, a significant reduction in the number of terms can be achieved 
by requiring that the LIV operators are rotationally symmetric. Let us then focus on 
the Lagrangian for neutrinos with LIV operators of mass dimension up to six involving 
a vector field u"- coupled to a Dirac neutrino ^/^ in a mass eigenstate with mass m 
that are quadratic in matter fields and hence modify the free field equations. We shall 
further assume that the Lorentz violation is diagonal in the mass basis. This might be 
justified by the idea that any theory of quantum gravity inducing such LIV must reduce 
to general relativity in the infrared, hence any Lorentz violation induced by quantum 
gravity would be primarily controlled by the charges that couple to gravity. Of course, 
this does not meant that the coefficients for each mass eigenstate are the same, as RG 
effects would not allow them to be the same at any energy. 

With these assumptions, the Lorentz violating terms (written in the mass basis) 
are exactly those for the QED fermions [1051 1120] 

— — 1 _ ^ 

- aiU^ij^Yi^i - biU^ij^'j^j'^i^i + -iCiU^Uyi)i'-i^ d" i)i (46) 

1 — ^ 1 — 



where i is the mass index. We've dropped the gauge covariant derivative above, as it 
is irrelevant and couples in the fiavor basis so merely would add needless complication. 
Also note that we have included both right and left projection operators, albeit, since 
standard model interactions only produce left-handed neutrinos, the constraints will 
primarily be on the corresponding left-handed operators. Finally, in contrast to the 
QED case, one cannot drop the Oj term here as this gives a contribution to the oscillation 
pattern. 

The above terms and the usual Dirac Lagrangian for the neutrino yield a high 
energy neutrino dispersion relation of 

E^ =p' + N^ (47) 

q (6) 4 

Nf = m\ + 2(a, + 6,)P - (q + d^V + 2(0,i - 0,2)^ + 2^^ 

As before, the term proportional to m/M-p\ has been ignored. The net result of the LIV 
terms is to modify the dispersion relation of each mass eigenstate according to 

^^-/-Kf = E^^^]i^' (48) 



n=l '"PI 

.(n) 



where ^^ is a coefficient that depends on the relevant terms in the Lagrangian (H711 . so 
that constraints on the ^^ can always be translated in constraints on the coefficients 
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of Eq. (H7|) . The corresponding anti-particle dispersion relation is easily derived by 
considering the behaviour of each term under CPT and given by 

E^-/-Kf = X:(-inr^^. (49) 

71=1 PI 

We leave the index n as a free phenomenological parameter and consider the cases 
n = 2,3,4 separately (the case n = 1 would produce huge effects at low energy and is 
strongly constrained). 

Note also that many existing neutrino oscillation experiments measure the 
transition probability between neutrino flavours Pu (/ and J here labelling neutrino 
flavours) which is affected by the modified dispersion relation fl47p . 

Pij = ^iJ-Yl 4^^^^. sin' ( -^ ] + 2Gij,, sin^ ( -^ j , (50) 

with SNfj = Nf — Nj and Fjjij and Guij are functions of the mixing matrixes. Many 
of these experiments also quote results on a deviation of the neutrino speed from that 
of light, i.e. 

/ A \ LIV 

(- = E~\6N^^ - 6ml) (51) 

which can be easily translated into constraints on the coefficients of Eq. (H7|) . 

5. Lorentz breaking and naturalness 

We presented in Section |3] an overview of several quantum gravity proposals/models 
which contributed in motivating the study of Lorentz breaking phenomenology. A 
common feature of these proposals is the presence of modified dispersion relations as 
those presented in the previous sections. For the rotational invariant case, these relations 
can be cast in the general form □ 

oo 

E^=p^ + m' + Y,V{n)P'' , (52) 

n=l 

where the fji are now dimensional coefficients. This is also why we considered increasing 
mass dimension operators in EFT providing dispersion relations at some definite order 
in n. 

However, while a definite QG model could predict different strengths for coefficients 
of different order (e.g. due too possible intrinsic symmetries of the model), from an EFT 
point of view the only relevant operators should be the lowest order ones, i.e. those of 
mass dimension 3,4 corresponding to terms of order p and p^ in the dispersion relation. 

+ Note that we are not considering here dissipative terms in the dispersion relation, as they would 
require a separate treatment (e.g. with respect to the preservation of unitarity) |121| . Nonetheless, 
such dissipative scenarios are logically consistent and even plausible within some quantum/emergent 
gravity frameworks. 
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Situations in which higher order operators "weight" as much as the lowest order ones 
are only possible at the cost of a severe fine tuning of the coefficients f/(„) . 

The reason for this is pretty simple: in EFT radiative corrections will generically 
allow the percolation of higher dimension Lorentz violating terms into the lower 
dimension terms due to the interactions of particles |122l 11231 I124J . Indeed, EFT loop 
integrals will be naturally cut-off at the EFT breaking scale, if such scale is as well the 
Lorentz breaking scale the two will effectively cancel leading to unsuppressed, coupling 
dependent, contributions to the mass dimension four kinetic terms that generate the 
usual propagators. Hence radiative corrections will not allow a dispersion relation with 
only p^ or p^ Lorentz breaking terms but will automatically induce extra unsuppressed 
LIV terms in p and p"^ which will be naturally dominant. One could argue that RG 
effects might naturally suppress the sizes of these coefficients at low energies. However, in 
specific models where the RG fiow has been calculated, the running of LIV coefficients is 
only logarithmic and so there is no indication that RG flow will actually drive coefficients 
to zero quickly in the infrared |116j . 

Current observational constraints are tremendous on dimension 3 operators and 
very severe on dimension 4 ones. This is kind of obvious, as dimension 3 operator would 
dominate at p — t- while the dimension 4 ones would generically induce a, species 
dependent, constant shift in the limit speed for elementary particles. Hence, within 
a "blind" EFT approach, it would deem unreasonable to further test departure from 
Lorentz invariance. 

Nonetheless, missing a definitive QG prediction the phenomenological approach 
that has been taken in the past is to fix an order of n and cast suitable constraints on 
the corresponding dispersion relation. This is is substantially equivalent to assuming a 
hierarchy of LIV coefficients of the sort 

^1 = ^1^, ^2 = r/2^, V3=V3j^ (53) 

or 

m = r/2^, V, = r/4^ (54) 

where M is the Lorentz breaking scale, that we assumed so far to be equal to the Planck 
scale, /i is some other far lower energy scale, and ?7(„) is now the de-demensionalized 
coefficient usually assumed to be 0(1). The exact hierarchy, whether it involves terms 
of every mass dimension as in ( l53l) . or only even dimension as in ( !54l) is model dependent. 
Several ideas have been advanced in order to justify such an unnatural (at least in 
EFT) hierarchy [82], sometime leading to a vibrant debate |125[I123] . but one can clearly 
see that the most straightforward solution for this problem would consist in breaking 
the degeneracy between the EFT scale and the Lorentz breaking one. We now briefiy 
describe threes ideas that have been put forward that would generate such scales. 
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5.1. A custodial symmetry? 

One (at least conceptually) simple solution, consists in appealing to some sort of 
"custodial symmetry", i.e. a symmetry, other than Lorentz, which would forbid the 
appearance of the lower dimension operators and is broken at the low energy (EFT) 
scale /i. The most natural symmetry that may play this role is of course supersymmetry 
(SUSY) [126[ I127J . Albeit SUSY is a symmetry relating fermions to bosons i.e. matter 
with interaction carriers, it also as a matter of fact, intimately related to Lorentz 
invariance. Indeed, it can be shown that the composition of at least two SUSY 
transformations induces space-time translations. However, SUSY can still be an exact 
symmetry even in presence of LIV and as such it can actually play the role of a custodial 
symmetry. 

It can be shown that the effect of SUSY on LIV EFT is to prevent dimension 
< 4, LIV operators to be present in the classical Lagrangian. Furthermore, the 
renormalization group equations for Supersymmetric QED with dimension 5 LIV 
operators as discussed in section 14.21 do not generate lower dimensional operators if 
SUSY is unbroken [M fT27j . 

Of course, SUSY is broken at low energies. The effect of soft SUSY breaking was 
also investigated in One should then expect, as it is confirmed by direct calculations 
[1261 1127] . that SUSY breaking allows the percolation of high energy LIV to the 
renormalizable operators. In particular, dimension k, ones arise from the percolation 
of dimension k + 2 LIV operators. (We consider here only k = 3, 4, for which these 
relationships have been demonstrated). However, in this case the coefficients appearing 
in front of the so generated low mass dimension operators will be naturally small 
being suppressed by powers of the ration of the SUSY breaking scale with the Planck 
scale. Typically the suppression is of order m^/Mpi (k = 3) or (ms/MpiY (k = 4), 
where rris is the scale of SUSY soft breaking. Although, given present constraints, the 
theory with k = 3 needs a lot of fine tuning to be viable, since the SUSY-breaking- 
induced suppression is not enough powerful to suppress enough linear modifications in 
the dispersion relation of electrons. Of course, one could always get rid e.g. in LIV QED 
of the terms of dimension 5 discussed in section l¥!2] by just imposing CPT invariance 
as a further symmetry of the problem and consider the next order operators like in 
section 14.31 . In this case one would end up inducing mass dimension 4 operators which 
are suppressed enough, provided m^ < 100 TeV. Current lower bounds from the Large 
Hadron Collider are at most around 950 GeV for the most simple models of SUSY [128] 
(the so called "constrained minimal supersymmetric standard model", CMSSM). 

As a final remark, let us notice that analogue models of gravity can be used as a 
particular implementation of the above mentioned mechanism for avoiding the so called 
naturalness problem via a custodial symmetry. This was indeed the case of multi- 
BEC [T29lfT30] . 
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5.2. Gravitational confinement? 

One possible alternative to the custodial symmetry idea consists in basically turning 
upside down the problem and assume that the Lorentz breaking scale is not set by the 
Planck scale. If one does this and begins with a theory which has higher order Lorentz 
violating operators only in the gravitational sector (e.g. Hofava-Lifshitz gravity), then 
one can hope that the gravitational coupling G^ ~ Mp^ will suppress the "percolation" 
to the matter sector where the constraints are strongest. Matter Lorentz violating terms 
will all possess factors of the order (/j,/Mpi)^ which can become strong suppression factors 
if /i <^ Mpi |131] . Of course, this is done at the dear price of large LIV effects in the 
gravity sector. However, the constraints we have on high energy violations in gravity 
are very poor so far, and hence this mechanism cannot be ruled out at the moment. We 
shall come back on this proposal later on in this review. 

5.3. Improved RG flows by strong dynamics 

One simple idea for solving the above mentioned naturalness problem is to construct 
a SME in such a way that the renormalization group flow of the LIV parameters (e.g. 
of 6c/c for elementary particles) will run to LI values in the IR. This was basically 
the original idea behind the pioneering work by Nielsen and collaborators in the early 
eighties [STJ |38l |39l SQl HH |12]. Unfortunately, as mentioned before, that program 
could never fulfil its scope because the LI values of the relevant quantities, albeit being 
attractors in the IR, they were so only with logarithmic fiow, i.e. too slow for preventing 
large deviations at relative low energies. However, there was recently a resurgence of 
activity in the field |132[ 1133] . The basic idea would be in this case to enhance at very 
high energies the RG flow via a strong coupled dynamics so that when approaching the 
logarithmic IR regime most of the suppression of LIV factors has already be achieved. 
Albeit preliminary, the results in this sense appear to be promising |133] . 

6. A toolkit for testing of Lorentz invariance 

Tests of Lorentz invariance have an early history given that tests of Special Relativity 
have been performed, in increasing accurate ways, sine its initial formulation in 1905. 
However, it is way more recent, as we discussed in our historical review, the development 
of a systematic EFT-based approach to the problem. While these kind of blind searchers 
are not directly related to a speciflc QG approach and hence can at time be difficult 
to asses in their physical consequences, it was undoubtedly an important step tho 
categorisation of the possible departures from extant Lorentz symmetry as they allowed 
a more clear understanding of what was going to be constrained and what kind of 
observations/experiments had to be seek for in order to do so. 

The phenomenological toolkit that was developed also thanks to this systematic 
approach is now quite rich and can be split in two big subsets: terrestrial experiments 
and astrophysical observations. 
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6.1. Experimental probes of low energy LIV: Earth based experiments 

It is an observational fact that Nature as we probe it well below the Planck scale it 
is Lorentz invariant to a very high degree. It is hence logic that in seeking for low 
energy deviations from this symmetry, as those systematically described by the minimal 
Standard Model extension (mSME), one has to resort to very high precision experiments. 
We shall now briefly list the most exploded observations in this sense (for more details 
see e.g. [TTIIITIITM]). 

6.1.1. Clock comparison Experiments Two co-local atomic transition frequencies can 
be considered as two clocks. As the clocks move in space, they pick out different 
components of the Lorentz violating tensors in the mSME. This in turn yields a sidereal 
drift between the two clocks which could be constrained by measuring the difference 
between the frequencies over long periods. This technique allows to cast very high 
precision limits on the parameters in the mSME (generally for for protons and neutrons). 

6.1.2. Cavity Experiments The technique adopted in cavity experiments casts 
constraints on the variation of the cavity resonance frequency (with respect to a 
stationary frequency standard) as its orientation changes in space. While this is 
intrinsically similar to clock comparison experiments, these kind of experiments allows to 
cast constraints also on the electromagnetic sector of the mSME (as one of the "clocks" 
in this case involves photons). 

6.1.3. Neutral mesons The mass difference of neutral mesons is one of the most 
accurately quantities in the SM. The SME operators do affect such a quantity in a 
Lorentz breaking way, hence one generically expects an orientation dependent change 
which can be constrained by looking for sidereal variations and other orientation effects. 
Also lifetime directional dependence has been investigated (see e.g. |135] ). 

6.1.4. Penning traps In a Penning trap a combination of static magnetic and electric 
fields confines a charged particle for long times. One can then problem possible 
deviation from exact Lorentz invariance by monitoring the particle cyclotron motion 
in the magnetic field and its Larmor precession due to the spin. In fact the relevant 
frequencies for both these motions are affected by some mSME operators and Penning 
traps can be set up so to make them very sensitive to differences in these frequencies. 

6.1.5. Spin polarized torsion balance Spin-torsion balances are a very effective tool for 
constraining the electron sector of the mSME. An example of such balances consists 
in an octagonal pattern of magnets which is constructed so to have an overall spin 
polarization in the octagons plane. Four of these octagons are suspended from a torsion 
fiber in a vacuum chamber so to give an estimated net spin polarization equivalent to 
K. 10^^ aligned electron spins. In order to detect Lorentz breaking effect one has again 
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to look for orientation dependent phenomena. For this reason the whole apparatus is 
mounted on a turntable. As the turntable moves, some of the Lorentz violation terms 
in the mSME, end up producing an interaction potential for non-relativistic electrons 
which induces a torque on the torsion balance. The torsion fiber is then twisted by an 
amount related to the relevant LIV coefficients which in this way can be constrained. 

The absence of a preferred direction is also checked with great precision using 
nuclear spin, which translates into more stringent limits on LIV operators for the 
light quarks and photons (the photon LIV operator will also contribute because of 
the electromagnetic interactions inside the nucleon). This possibility was analysed for 
the mSME as week as for the SME with dimension 5, CPT odd operators [85]. These 
experiments are again based on a search for orientation dependent phenomena using 
comagnetometers of different types [136t I137J able to detect the sidereal variation with 
great accuracy. These methods are currently able to constraints the coefficients quarks 
and photons at 0(10"^ ^ IQ-^). 

Remarkably, the constraints of O(10~^^) cast in |137j on anisotropic components of 
the limit speed for neutrons in the mSME can also be converted in a severe constraint 
also on any time-like component (i.e. rotational invariant). In fact, we assumed the 
preferred frame to be aligned with the CMB but the solar system is moving relative to 
such rest frame with boost factor ~ 10^'^ at a declination of —7°. The missing alignment 
will hence always induce an anisotropic component suppressed by the square of the 
boost factor. In conclusion, this allows to constrain an eventual isotropic component in 
the neutron limit speed at 0(10"^"^). 

6.1.6. Recoil of slow atoms In alternative to the relativistic expansion that we have 
considered so far, one might as well consider non-relativistic limit, p -^ m, oi the 
phenomenological dispersion relations of the form (l52l) supplemented with the further 
assumption fl53|) . It is easy to see that in such a case the dispersion relations and its 
most relevant correction would look like 

V 1 / v \ 

E ^ m + h 7777 Vi'iTT'P + V2P'^ H • (55) 



2m 2M y ' m^ 

It was recently noticed |138[ 1139] that ultra-precise cold-atom-recoil experiments have 

now reached sufficient sensitivities for constraining the lowest coefficients in the above 

expansion (albeit not strongly). The basic techniques consists in using the available 

measurements of the "recoil frequency" of atoms with experimental setups involving 

one or more "two-photon Raman transitions" |140lll411ll42] . The basic idea is that one 

can give momentum to an atom through a process involving absorption of a photon of 

frequency v and a stimulated emission, in the opposite direction, of a photon of frequency 

v' . A straightforward application of energy- momentum conservation implies |138[ 1139] 

/\v h 
= (5g) 

2i^*(z/* +p/h) m 
where i^^, is resonance frequency of the atom, p is the momentum transferred to the 
atom and Av ^ v — v' . The crucial point is that hlS.v ~ E(j) + hu + hv') — E{p) ~ 



Tests of Lorentz invariance 24 

E{2hu^+p)—E{p) and hence would be affected by modification of the dispersion relation, 
modifications which, thanks to the high precisions of current measurements could be 
then use to constraint the lowest order coefficients rji and 772 [13811139] . 

6.2. Observational prohes of high energy LIV: astrophysical QED reactions 

Testing the higher mass dimensions operators (5 or 6) of the SME is obviously a task that 
requires to probe much higher energies than those achieved in terrestrial experiments. 
Energies up to 10^*^ eV are achieved in high energy astrophysics and as such this field 
has played an eminent role in QG phenomenology. One comment is in order before we 
further describe these tests and it concerns our parametrisation. We have introduced 
dimensionless coefficients t]*-"' and expressed the higher order terms in the dispersion 
relations via suitable rations of the particle momentum to the Planck scale. Of course 
nothing guarantees that models of quantum/emergent spacetime will predict the Loretz 
breaking scale to be coincident with the Planck scale. So, missing a derivation of our 
dispersion relation from a specific QG model, this has to be considered just a convenient 
parametrisation. 

Let us begin with a brief review of the most common types of reaction exploited in 
order to give constraints on the QED sector. For definiteness, we refer to the following 
modified dispersion relations: 

E',=e + &^^, Photon (57) 
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E^i = m^ + p + 77± — ^^32 Electron-Positron , (5J 

iVj[r)l 



where, in the EFT case, we have ^^"^^ =(,).= (— )"C- ^^'i V = V+ (— )"^- 

6.2.1. Photon time of flight Photon time-of-flight were historically among the first 
tests proposed |19] for QG phenomenology and although within an EFT framework 
they currently provide limits several orders of magnitude weaker than the best ones, they 
have been widely adopted in the astrophysical community. Furthermore, being purely 
kinematical in nature, they apply on any modified dispersion relation independently 
from the underlying dynamical framework. Let us then discuss a general description 
of time-of-flight effects, elaborating later about their application to the specific case of 
EFT. 

A photon dispersion relation in the form of ( 1571) implies that photons of different 
colours (wave vectors ki and ^2) travel at slightly different speeds. We shall assume 
here that no birefringent effects are present, so that ^_^ = ^_ (or alternatively in EFT 
we shall consider one helicity at a time). We propagating on a cosmological distance d, 
the effect of energy dependence of the photon group velocity will generically produces 
an energy dependent time delay 

1 i,'^"^ l,'^~2 

AtW = ^ "' ~ 7 e^") d , (59) 
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which increases hnearly with dij and with the energy difference as long as n > 2. The 
largest systematic error affecting this method is the uncertainty about whether photons 
of different energy are produced simultaneously in the source (generally Gamma-Ray 
Bursts or Active Galactic Nuclei). 

So far, the most robust constraints on ^^'^\ derived from time of flight differences, 
have been cast adopting a statistical analysis applied to the arrival times of sharp 
features in the intensity at different energies from a large sample of GRBs with known 
redshifts [1431 1144] . leading to limits ^^^^ < O(IO^). A recent example illustrating the 
importance of systematic uncertainties can be found in [145j . where the strongest limit 
^(3) ^ 47 jg found by looking at a very strong flare in the TeV band of the AGN Markarian 
501. Furthermore, analysing the flare time-structure, the authors also obtained a best 
fit for $^^^' ~ 0(1). However, a similar fit could also be achieved by standard plasma 
physics. 

In order to alleviate systematic uncertainties one could resort to extra information 
based on the specific dynamical model adopted. In particular in EFT one expects 
dimension 5 CPT odd terms to produce birefringent photon dispersion relation. One 
could then try to measure the velocity difference between the two polarization states at 
a single energy photon, corresponding to 

At = 2\^^^^\kd/Mpi. (61) 

This bound would require that both polarizations be observed and that no spurious 
helicity-dependent mechanism (e.g. propagation through a birefringent medium) affects 
the relative propagation of the two polarization states. As such this possibility has not 
been exploited so far. 

More interestingly, Eq. (l59l) is not necessarily valid in birefringent theories. In fact, 
photon beams generally are not circularly polarized; thus, they are a superposition of 
fast and slow modes. The net effect of this superposition may end-up erasing the time- 
delay effect. We can compute this effect on a generic photon beam in a birefringent 
theory by means of the associated electric field, and let us assume that this beam has 
been generated with a Gaussian width 

E = A L^i^ot-k+{no)z) ^-{z-v+t)^5nlg^ _^ ^iiUot-k- ino)z) ^-(z-v-t)^5nlg\ ^ ^q2\ 

where Qq is the wave frequency, SQq is the gaussian width of the wave, k'^^Qo) is 
the "momentum" corresponding to the given frequency according to (1571) and e± = 
(ei lb ie2)/y/2 are the helicity eigenstates. Note that by complex conjugation e^ = e_. 
Also, note that k'^{u) = w =F ^u'^/Mpi. Thus, 



* 



Note that for an object located at cosmological distance (let z be its redshift), the distance d becomes 

1 /■" 1+2' 
d{z) = — / dz' , (60) 

where d{z) is not exactly the distance of the object as it includes a (1 + z)^ factor in the integrand to 
take into account the redshift acting on the photon energies. 
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The intensity of the wave beam can be computed as 

E ■ E* = Up ("e^^^^o/A^pi^ + ^-2iinl/Mpiz\ ^-5nl((z-v+tf+(z-v-tf) 

= 2|Ape-2^^o(^-*)^ cos (2^^noz] e'^^'^^^^"*)' . (64) 

We see that there is a modulation of the wave intensity that depends quadratically on 
the energy and hnearly on the distance of propagation. In addition, for a Gaussian wave 
packet, there is a shift of the packet centre, that is controlled by the square of ^^^'^ /Mpi 
and hence is strongly suppressed with respect to the cosinusoidal modulation. So far, 
to the best of our knowledge, this EFT specific features have not been tested (and as 
we shall see below they would be extremely difficult to see anyway as birefringence is 
strongly constrained by current observations). 

6.2.2. Vacuum Birefringence The fact that in the "LIV extended QED" discussed in 
section 14.21 opposite "helicities" have slightly different group velocities, implies that the 
polarisation vector of a linearly polarised plane wave with energy k rotates, during the 
wave propagation over a distance d, through the angle 

^(d) = ""^^^:""^^^ ci^e^3) fc!^. (65) 

^ ^ 2 ^ 2Mpi ^ ^ 

It is obvious that such angle will be different for different photon energies hence this effect 

could potentially disrupt the amount of polarisation present in a some polarised light 

travelling over long distances. More specifically, depending on the amount of available 

information on the specific object observed, one can use this fact to cast constraints on 

the photon LIV parameter |^*^^)| in two different ways |146] : 

Decrease in polarization degree. If some net amount of polarization is measured in 
the band say ki < E < k2, a.n order-of- magnitude constraint arises from the fact that 
the angle of polarization rotation f l65|) cannot differ by more than 7r/2 over this band, 
as in this case the detected polarization would fluctuate sufficiently for the net signal 
polarization to be suppressed |147l 1117] . From ( l65l) . assuming A^ < 7r/2 implies 



vrMpi 
m - kl)d{z) 



e^ < n/,L.. , (66) 



This constraint merely relies on the detection of a polarized signal and does not make 
use of the observed polarization degree but one can obtain for a more refined limit by 
calculating the maximum observable polarization degree, given some maximum intrinsic 
value [TI8] : 

n(0 = n(0) y/(cos(2^))2, + (sin(2^))2„ (67) 

where 11(0) is the maximum initial degree of polarization, 6 is defined in Eq. fl65|) and the 
average is weighted over the source spectrum and instrumental efficiency, represented 
by the normalized weight function V{k) |147] . Conservatively, one can a priori set 
n(0) = 100%, however for many astrophysical sources a lower value may be justified 
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on the basis of theoretical inodelhng. Using (1^7|) . one can then cast a constraint by 
requiring n(,^) to exceed the observed value. 



Rotation of polarization angle. A different possibility for casting a constraint consists 
in looking at the overall rotation of the average polarisation. Let's suppose that one 
observes polarised light over some energy band which shows an overall position angle 
^obs with respect to a fixed direction. At fixed energy, the polarization vector rotates by 
the angle f l65|) (Faraday rotation is negligible at high energies.); if the position angle is 
measured by averaging over a certain energy range, the final net rotation (A^) is given 
by the superposition of the polarization vectors of all the photons in that range: 

(sin(20))^ 



tan(2(A^)) 



(6J 



(cos(2^))^ ' 

where 6 is given by fl65|l . If the position angle at emission 6^ in the same energy band is 
known from a model of the emitting source, a constraint can be set by imposing 

tan(2 (A^)) < tan(2^obs - 2^;) . (69) 

While this limit is tighter than those based on eqs. (166|) and ( 167|) . it clearly hinges 
on assumptions about the nature of the source, which may introduce significant 
uncertainties. 

In conclusion the fact that polarised photon beams are indeed observed from distant 
objects imposes strong constraints on LIV in the photon sector (i.e. on ^^^^), as we shall 
see later on. 



6.2.3. Threshold reactions LIV corrections are quite important in threshold processes 
because the LIV term (which as a first approximation can be considered as an additional 
mass term) only needs to be comparable to the (invariant) mass of the particles produced 
in the final state for strongly affecting this kind of reactions. I.e. a rough estimate for 
the threshold energy is 



Pth 



rri^M, 



n~2\ 1/" 



PI 
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(70) 



where m is the typical mass of particles involved in the reaction. Interesting values 
for pth are discussed, e.g., in [20] and given in Tab. [U Reactions involving neutrinos 



Table 1. Values of pth, according to cq. ([70)1 . for different particles involved in the 
reaction: neutrinos, electrons and proton. Here we assume 77^"^ ~ 1. 





nil, — 0.1 eV 


rUe ~ 0.5 MeV 


mp ~ 1 GeV 


n = 2 


0.1 eV 


0.5 MeV 


1 GeV 


n = 3 


500 MeV 


14TeV 


2 PeV 


n = A 


33TeV 


74PeV 


3EeV 
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are the best candidate for observation of LIV effects, whereas electrons and positrons 
can provide results for n = 3 theories but can hardly be accelerated by astrophysical 
objects up to the required energy for n = 4. In this case reactions of protons can be 
very effective, because cosmic-rays can have energies well above 3 EeV. 

For the above reason a quite rich phenomenology of threshold reactions is introduced 
by LIV in EFT and threshold theorems can be generalized [81]. Sticking to the 
present case of rotational invariance and monotonic dispersion relations (see |149] for 
a generalization to more complex situations), the main conclusions of the investigation 
into threshold reactions are that |80] 



• Threshold configurations still corresponds to head-on incoming particles and 
parallel outgoing ones 

• The threshold energy of existing threshold reactions can shift, and upper thresholds 
(i.e. maximal incoming momenta at which the reaction can happen in any 
configuration) can appear 

• Pair production can occur with unequal outgoing momenta 

• New, normally forbidden reactions can be viable 

Let us now briefly review the main reaction used so far in order to casts constraints. 

6.2.4- LIV-allowed threshold reactions: '-/-decay While normally the decay of a photon 
in a lepton-antilepton pair is kinematically forbidden, it might be allowed in the presence 
of LIV induce modified dispersion relations of the sort we considered so far. Taking the 
case n = 3 and ^'^^•* = as a simple example, the threshold is set by the condition [82] 

A:,, = (6v^m^Mpi/|r/l^)|)i/\ (71) 

Noticeably, as already mentioned above, the electron-positron pair can now be created 
with slightly different outgoing momenta (asymmetric pair production) [801 181] . 
Furthermore, the decay rate is extremely fast above threshold [H2| and is of the order 
of (10 ns)-^ {n = 3) or (10"^ ns)^^ (n = 4). 

6.2.5. LIV-allowed threshold reactions: Vacuum Cerenkov and Helicity Decay In the 
presence of LIV, also the process of Vacuum Cerenkov (VC) radiation e^ — ?■ e^7 
can occur (this can be easily seen, as for the gamma decay, by imposing the energy 
momentum conservation). If we set again n = 3 and ^^^'^ ^ and , the threshold energy 
is given by 

pvc = (m^Mpi/2r/(3))^/^ ~ 11 TeV r/^^/^ . (72) 

Just above threshold this process is also extremely efficient, with a time scale of order 
rvc ~ 10-9 s [12] . 

A slightly different version of this process is the Helicity Decay (HD, e=p — )■ e±'-f). 
If rj^ ^ f]-, an electron/positron can flip its helicity by emitting a suitably polarized 
photon. This reaction does not have a real threshold, but rather an effective one [82] 
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— Pud = {mlMpi/ ArjY'^ , where Ar] = \r]\, — ''7_ | — at which the decay hfetime thd 
is minimized. For Arj ^ 0(1) this effective threshold is around 10 TeV. Note that 
below threshold Thd > A?7~^(p/10 TeV)~^10~^s, while above threshold Thd becomes 
independent of Arj p2] . 

Apart from the above mentioned examples, one can also look for modifications of 
normally allowed threshold reactions especially relevant in high energy astrophysics. 

6.2.6. LIV-allowed threshold reactions: photon splitting and lepton pair production 
Once photon decay and vacuum Cerenkov are allowed also the related relations in which 
respectively the our going lepton pair is replaced by two or more photons, 7 — )■ 27 
and 7 — 7> 37, etc. , or the outgoing photons is replaced by an electron-positron pair, 
e~ — )■ e^e^e"*", are also allowed. 

Photon splitting It is easy to show that this process requires UV superluminal photons, 
i.e. ^'^"•* > [80]. When allowed, the relevance of this process is simply related to its 
rate. The most relevant cases are 7 — ;■ 77 and 7 — )■ 87, because processes with more 
photons in the final state would be suppressed by more powers of the fine structure 
constant. 

The 7 — ;■ 77 process is forbidden in standard QED because of kinematics and 
CP conservation. In LIV QED neither condition holds, however, we can argue that 
this process is suppressed by an additional power of the Planck mass, with respect to 
7 — )■ 87. In fact, in LI QED the matrix element is zero due to the exact cancellation 
of fermionic and anti-fermionic loops. In LIV EFT this cancellation is not exact and 
the matrix element is expected to be proportional to at least ^(E'/Mpi)^, p > 0, as it is 
induced by LIV and must vanish in the limit Mpi — ?> 00. 

Let us consider then 7 — )■ 87. This process has been studied for n = 8 in 
[8OI 1150] . In particular, in |150] it was found that, if the "effective photon mass" 
m'i = C,E!^/Mp^^ <^ ml, then the splitting lifetime of a photon is approximately 
^n=3 ^ 0.025 ^~^/~^ (50 TeV/E^) s, where / is a phase space factor of order 1. This 
rate was rather higher than the one obtained via dimensional analysis in [80] because, 
due to integration of loop factors, additional dimensionless contributions proportional 
to TTig enhance the splitting rate at low energy. 

Lepton pair production The process e~ — )■ e~e~e+ is similar to vacuum Cerenkov 
radiation or helicity decay, with the final photon replaced by an electron-positron pair. 
While various combinations of helicities for the different fermions can be considered, we 
shall consider here, for illustrative reasons, the most simple case where all electrons have 
the same helicity and the positron has the opposite helicity, so that the threshold energy 
will depend on only one LIV parameter. The threshold for this reaction was derived in 
In [80] and found higher than that for soft vacuum Cerenkov radiation by a factor 2.5. 
Given that the reaction rate is high as well, one can casts constraints by using just the 
value of the threshold. 
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6.2.7. LlV-modified threshold reactions: Photon pair-creation A process related to 
photon decay is photon absorption, 770 — > e~^e~. Unhke photon decay, this is allowed 
in Lorentz invariant QED and it plays a crucial role in making our universe opaque to 
gamma rays above tents of TeVs as these are absorbed in this way by the infra-red and 
cosmic microwave background radiation. 

Let's call cuq the (low) energy associated to the target photon 70, then the threshold 
for the reaction occurs in a head-on collision with the second photon having momentum 
khi = ni? /ujq. For example, if /cli = 10 TeV (the typical energy of inverse Compton 
generated photons in some active galactic nuclei) the soft photon threshold Wq is 
approximately 25 meV, corresponding to a wavelength of 50 microns. 

In the presence of Lorentz violating dispersion relations the threshold for this 
process is in general altered, and the process can even be forbidden. Moreover, as 
we mentioned before, in some cases an upper threshold con develop for which the 
process does not occur |15H [80 | |8T |I149] . Physically, this means that at sufficiently high 
momentum the photon does not carry enough energy to create a pair and simultaneously 
conserve energy and momentum. Note also, that an upper threshold can only be found 
in regions of the parameter space in which the 7-decay is forbidden, because if a single 
photon is able to create a pair, then a fortiori two interacting photons will do [80] . 

Let us exploit the above mentioned relation 77^ = (— )"'?7t between the electron- 
positron coefficients, and assume that on average the initial state is unpolarized. In this 
case, using the energy-momentum conservation, the kinematics equation governing pair 
production is the following [82] 

i?^> - ^ + s^ - sr <^^"'" + '-'" '^ - ^'" ' <''' 

where < y < 1 is the fraction of momentum carried by either the electron or the 
positron with respect to the momentum k of the incoming high-energy photon. In 
general the analysis is rather complicated. In particular it is necessary to sort out 
whether the thresholds are lower or upper ones, and whether they occur with the same 
or different pair momenta 

6.2.8. Synchrotron radiation Synchrotron emission is strongly affected by LIV [152[ 
|52| I153[ I154[ |82| I155[ I156[ 1157] , however for Planck scale LIV and observed energies, it 
is a relevant "window of opportunity" for QG phenomenology only for dimension four 
or five LIV QED albeit we shall see that in specific QG models (see below section |9^ . 
with the Lorentz breaking scale not equal a priori to the Planck one, it might play a 
relevant role also for dimension six LIV operators. For this reason we shall generically 
indicate the Lorentz breaking scale Muv in what follows. 

In both LI and LIV cases [82], most of the radiation from an electron of energy E 
is emitted at a critical frequency 

^. ^ leB^ (74) 
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where 'j{E) = (1 — v'^{E))~^/'^, and v{E) is the electron group velocity. Assuming 
Hamiltonian dynamics, the group velocity of an electron with a modified dispersion 
relation at arbitrary order as in ( 1571) can be computed as (neglecting higher order 
terms) |157] 

We immediately see that v{E) can exceed 1 if r/ > or it can be strictly less than 1 if 
?7 < once the LIV term is prevailing on the suppressed mass term. This introduces a 
fundamental difference between particles with positive or negative LIV coefficient rj and 
the possibility to see the effect even for p/Muy -C 1. In [152j a simple, but stringent, 
constraint for n = 3 was cast using the fact that subluminal dispersion relations (77 < 0) 
admit a maximal critical frequency for the synchrotron emission. For arbitrary n > 2 
this would be |157] 

2m V 3n) \-r]{n-l){3n-4:)) ' ^ ' 

Basically, if synchrotron emission up to some frequency Wobs is observed, one can deduce 
that the LIV coefficient for the corresponding leptons cannot be more negative than the 
value for which uj"^'^^ = oJobs- 

If T] is instead positive the leptons can be superluminal. One can show, e.g. for 
n = 3 and Mliv = MpL, that at energies Ec^8 TeV/rj^^^, liE) begins to increase faster 
than E/rUe and reaches infinity at a finite energy, which corresponds to the threshold for 
soft VC emission. The critical frequency is thus larger than the LI one and the spectrum 
shows a characteristic bump due to the enhanced Uc- While this is the general behaviour 
to be expected, a derivation of a constraint in this case requires a direct confrontation 
with observation and hence detailed reconstruction of the synchrotron spectrum in the 
presence of LIV effects |156[ 1157] . 

7. Experimental and observational constraints 

We shall now give a brief overview of the most stringent constraints available on the 
different sector of the SME. Of course, different operators, of different mass dimension, 
require different sort of observations and hence it will make sense to divide our discussion 
not only by SM sector but also by the order of the operators considered. 

7.1. Rotational invariant QED in mSME 

Lorentz and CPT breaking operators does not depend on energy in the mSME. Therefore 
it is generally true that high precision laboratory and terrestrial experiments might be 
preferable with respect to astrophysical observation in this case. For a recent listing of 
all the constraints including the non-rotationally invariant case, see |134] , we shall limit 
ourself here to the rotationally invariant case discussed before. 
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7.1.1. Spin polarized torsion constraints on electrons: Spin polarized torsion balances 
can place limits on the electron sector of the niSME [158] and can be very roughly 
described as a magnet attached on a torsion fiber. Hence, torsion balances are 
characterised by a large number of aligned electron spins. The mSME coefficients 
give rise to an interaction potential for non-relativistic electrons which produces an 
orientation dependent torque on the torsion balance which can be measured using the 
twist of the torsion fiber. The torsion balance in its sealed vacuum chamber is mounted 
on a rotating turntable, which allows for very sensitive detection of any anomalous 
torque as a function of the rotation frequency (and the Earth's rotation and motion 
in the solar system). These experiments sets limits on the b coefficient in the electron 
modified Lagrangian Eq. (I32l) . Current limits are of the order of 10~^^ GeV [158] . 

7.1.2. Accelerator constraints on electrons: Energy loss mechanisms such as the 
vacuum Cerenkov effect g — ?■ g + 7 (see section I6.2.3P for charged fermions are 
incompatible with a constant beam energy as long as the energy loss rate is high enough. 
The LEP experiment accelerated electrons and positrons to energies of roughly 100 GeV 
in the lab frame and it is known that the beam does not lose any significant energy to 
vacuum Cerenkov radiation. Additionally the synchrotron emission from the beam has 
been measured. Using the characteristics of the LEP beam and the related synchrotron 
radiation spectrum sets limits on the coefficient c for electrons (as in Eq. ( 132|) ) to be of 
order \c\ < 10~^^. 



7.1.3. Astrophysical inverse Compton bounds on electrons: The parameter din Eq. (13211 
controls a spin dependent change to the dispersion relation for a fermion. We observe 
energies from the radio up to 80 TeV from severe active galactic nuclei as well as pulsar 
wind nebula, e.g. the Crab nebula [159] . The overall spectrum of these sources is well 
understood and the high energy emission is dominated by inverse Compton scattering of 
accelerated electrons off of lower energy photons (mainly synchrotron and CMB photons 
depending on the specific astrophysical object). The magnetic field present in the same 
sources causes the electrons to precess, thereby destroying any initial polarization of the 
electrons. Therefore, one can argue that both helicities of electrons must be present and 
stable if there is to be inverse Compton radiation |160] . This fact forbids the vacuum 
Cerenkov effect for both helicities, implying a double sided constraints on d of the order 
of |rf| < 0(10-^2) |X6o] in the CMB rest frame. 

7.1.4. Cosmic ray and HESS bounds on photons: The observation of high energy 
cosmic ray protons and TeV gamma rays allows to cast constraints on the on kp 
parameter for photons in Eq. ( 133|) . Assuming unmodified protons (within the mSME 
which sector is the "unmodified" sector is arbitrary, in that the limiting speed of one 
of the sectors can be defined to be the "speed of light"), in [161] the authors used 
the necessary absence of vacuum Cerenkov p — )■ p + 7 and a Pierre Auger event in 
conjunction with the excess of high energy TeV gamma rays observed by the HESS 
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telescope (which forbids the gamma decay 7 — )■ p + p) to produce a two sided bound on 
kp of -9 X 10-16 < A;^. < 6 X IQ-^o. 

7.1.5. Colemann-Glashow limit of the mSME QED: One can of course restrict the 
QED sector in the mSME to the rotational invariant subset and solely dimension four 
operators. In this case the model basically coincides with the Coleman-Glashow one |55j . 
Also in this case the constraints are quite strong, for example on the QED sector one can 
easily see that the absence of gamma decay up to 50 TeV provides a constraint of order 
IQ-i^ on the difference between the limit speed of photons and electrons |119j . Similarly, 
the observations of protons up to 10^^ eV in cosmic rays is sufficient evidence for a 
constraint of order 10"^° on the difference of speed between them and light (constraints 
of order 10"^^ -=- lO"^'^ can also be achieved, and are sometimes quoted in the literature, 
if one takes slightly higher values for the observed energies of propagating protons. We 
report here a more conservative constraint given the present uncertainties on the highest 
energy tail of the cosmic rays spectrum, see Section [775|) . 

7.2. QED with rotational invariant CPT odd dimension 5 ops 

There is a very rich literature devoted to constraining the Myers-Pospelov version of the 
SME. We shall not attempt here a detailed discussion but just briefly summarise the 
available constraints. 

7.2.1. Synchrotron radiation constraints from the Crab Nebula: Presently a very 
important object in this sense is the Crab nebula (CN). This pulsar wind nebula had 
its origin by a supernova explosion observed in 1054 A.D. Its distance from Earth is 
approximately 1.9 kpc and it is one of the best studied sources of diffuse radio, optical 
and X-ray radiation. The Nebula emits an extremely broad-band spectrum (21 decades 
in frequency, see |156] for a comprehensive list of relevant observations) that is produced 
by two major radiation mechanisms. The emission from radio to low energy 7-rays 
{E < 1 GeV) is thought to be synchrotron radiation from relativistic electrons, whereas 
inverse Compton (IC) scattering by these electrons is the favoured explanation for the 
higher energy 7-rays. From a theoretical point of view, the current understanding of 
the whole environment is based on the so called Kennel-Coroniti |162] . which accounts 
quite accurately for the general features observed in the CN spectrum. 

Developing on the initial idea discussed in |152] . a full reconstruction of the 
synchrotron emission processes in the CN within the SME with dimension 5 CPT 
operators (see Eqs. (1371) and ( l38l) ) has been performed in |156] . This procedure requires 
fixing most of the model parameters using radio to soft X-rays observations (which are 
basically unaffected by LIV), a careful step by step revisitation of the basic assumptions 
made in the Kennel-Coroniti model and also to take into account the contribution of 
novel LIV effects such as vacuum Cerenkov and helicity decay. Finally, a x^ analysis 
can be performed to quantify the agreement between models and data |156] . From this 



Tests of Lorentz invariance 34 

analysis, one can conclude that the LIV parameters for the leptons are both constrained, 
at 95% CL, tobe |r/±| < 10"^ 

7.2.2. Birefringence constraints: The best constraints on the photon sector are instead 
obtained by using birefringence effects associated with the CPT odd nature of the 
relevant LIV operators. Strong constraints came again from the CN |146j . from which a 
value l^^'^^l < 6 X 10~^° at 95% Confidence Level (CL) was obtained by considering the 
observed polarization of hard-X rays |163] (see also |164] ). Polarized light from GRBs has 
also been detected and given their cosmological distribution they could be ideal sources 
for improving the above mentioned constraints from birefringence. Attempts in this 
sense were done in the past |117l I165J but the observed of polarisation was later deemed 
controversial. Furthermore, GRBs for which the polarization is detected and the spectral 
redshift is precisely determined are scarce. In |166] this problem was circumvented 
by using indirect methods (the same used to use GRBs as standard candles) for the 
estimate of the redshift. This leads to a possibly less robust but striking constraint 
1^(3) I < 2.4 X lO""*^^. A similar constraint |.^*^'^^| < 10"""^^ was obtained in [167j by making 
use of two GRBs observed by the Japanese satellite IKAROS. Also in this case the 
redshift was determined using indirect methods such as the correlation between peak 
energy and luminosity of the GRB prompt emission. Remarkably, the above mentioned 
constraints were recently improved by using the INTEGRAL/IBIS observation of the 
GRB 061122, for which a redshift z ~ 1.33 was derived thanks to the determination of 
the GRB's host galaxy. In this case a constraint \^^'^'^\ ^ 3.4 x 10~^^ was derived [168] 
(see also |169j for a similar previous constraint). 

7.3. QED with rotational invariant CPT even dimension 5 and 6 ops 

As we have shown before, this kind of operators lead to modified dispersion relations 
characterised by quartic terms in the particle momenta suppressed by the squared Planck 
mass. Furthermore the CPT even nature of the operators prevent birefringent features 
in the photon sector. This feature basically implies the impossibility to use any of 
the above discussed observation for casting effective constraints (e.g. the synchrotron 
emission from the CN would not provide a constraint better that 10^ on the electron 
parameter given the [E/M-pif suppression). However, there are available observations 
exploring much higher energies, i.e. those related to the so called Ultra High Energy 
Cosmic Rays (UHECR). 

7.3.1. Astrophysical constraints from GZK reaction secondaries: The GZK cut off 
[98| [99] , is a suppression of the high-energy tail of the UHECR spectrum arising from 
interactions with CMB photons: p7 — ?> A'^ — t- p-K^irm^). When Lorentz invariance 
holds, this process has a threshold energy Ei^a — 5 x 10^^ (ci;o/1.3 meV)~^ eV (where 
Uq is the target photon energy). Experimentally, the presence of a suppression of the 
UHECR flux was claimed only recently [SI Il3] (and, as discussed in section [3l it was 
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initially claimed to be absent by the AGASA collaboration |1UU] ). Although the cut off 
could be also due to the finite acceleration power of the UHECR sources, the fact that 
it occurs at the expected energy favours the GZK explanation. The results presented 
in |170j seemed to further strengthen this hypothesis (but see further discussion in the 
conclusions) . 

Significant limits on ^ = ^^^^ and r] = r]^^^ for the electron/positron can be 
derived by considering UHE photons generated as secondary products of the GZK 
reaction |17H 1172] . These UHE photons originate because the GZK process leads to the 
production of neutral pions that subsequently decay into photon pairs. These photons 
are mainly absorbed by pair production onto the CMB and radio background. Thus, 
the fraction of UHE photons in UHECRs is theoretically predicted to be less than 1% 
at 10^^ eV |173j . Several experiments imposed limits on the presence of photons in the 
UHECR spectrum. In particular, the photon fraction is less than 2.0%, 5.1%, 31% and 
36% (95% C.L) ai E = 10, 20, 40, 100 EeV respectively [Hll [175] (at lower energies 
even stricter limits 0.4%, 0.5%, 1.0%, 2.6% and 8.9% (95% C.L) above E = 1, 2, 3, 5 
and 10 EeV respectively were recently found |176] ). 

The key idea for casting constraints here, relies on the strong dependence of the 
photon pair production on on LIV modifications (see section I6.2.7p . In particular, 
the (lower) threshold energy can be shifted and in general an upper threshold can be 
introduced [80] . If the upper threshold energy is lower than 10^° eV, then GZK secondary 
UHE photons are no longer attenuated by the CMB and can reach the Earth. Hence, 
they would constitute a significant fraction of the total UHECR flux and thereby they 
would violate the above mentioned experimental bounds |17H I172[ 1177] . Of course, 
one LIV is introduced these extra UHE photons could now be removed via the above 
discussed photon splitting process. However, the analysis of the rate we presented in 
|150] and discussed in section 16.2. 6[ does not apply to photons around 10^^ eV given that 
at these energies m? ^ ml ii C,^^^ > 10"^'^ and ^^^^ > 10"^. In this regime, contributions 
to the rate are in this case at most logarithmic, as the momentum circulating in the 
fermionic loop is much larger than rrie- Moreover, the splitting rate depends only on m^, 
the only energy scale present in the problem. One then expects the analysis proposed 
in [HO] to be correct and the splitting time scale to be negligible at E^ ^ 10^^ eV. 

To complete the analysis one can also notice that the 7-decay process can also imply 
a significant constraint. Indeed, if some UHE photon {E^ ~ 10^^ eV) is detected by 
experiments (and the Pierre Auger Observatory, PAO, will be able to do so in few years 
|174j ). then 7-decay must be forbidden above 10^^ eV |172] (note that for n = 4 the two 
photons felicities travel at the same speed, so there is in this case a single population). 

The combination of above mentioned "upper threshold constraint" and the expected 
constraint from the absence of gamma decay would bound the QED LIV parameters at 
order n = 4 to the roughly rectangular region (see Figure 2, of |178] for more details) 

- 10"' < e^^^ < 10"^ - 10"' < r/W < 10-6 (77) 
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7.4- SME: Hadronic sector 

The very same GZK photo-pion production is strongly affected by LIV and the 
constraints imphed by the detection of this effect have been extensively considered in 
the literature [1791 [1801 EOl W, [THU [I82]. Nonetheless, a detailed LIV study of the GZK 
feature is hard to perform, because of the many astrophysical uncertainties related to 
the modelling of the propagation and the interactions of UHECRs. In fact, the Lorentz 
breaking operators could lead to relevant modifications of the GZK mean free path. 
Consequently, the propagated UHECR spectrum can display new features, hke bumps at 
specific energies, suppression at low energy, recovery at energies above the cutoff. These 
are all features which cannot be easily conciliated with the observed spectrum, even 
taking into account experimental uncertainties. Furthermore, the emission of Cherenkov 
7-rays and pions in vacuum would lead to sharp suppression of the spectrum above the 
relevant threshold energy. After a detailed statistical analysis of the agreement between 
the observed UHECR spectrum and the theoretically predicted one in the presence of 
LIV and assuming pure proton composition, the final constraints implied by UHECR 
physics are (at 99% CL) [183] 

- 10"' < # < 10-^ 

-10-' < V^:^ < 10-^ (r/W > 0) or < 10"^ (r/^) < 0) . (78) 

Obviously the same analysis can be applied also to dimension five, CPT odd, operators 
leading to much stronger constraints (order O(10~^'^)). 

7.5. An open issue from UHECR 

We just saw how the GZK cutoff in UHECR physics has played a determinant role 
in casting constraints for n = 4 (CPT even dimension 5 and 6 operators) dispersion 
relations in the QED and Hadronic sectors of the SME. It is hence quite unfortunate 
that these constraints are still suffering from experimental uncertainties about the 
composition of the UHECR and hence the nature of the steep fall-off in the spectrum 
observed above 10^^ eV. Indeed, UHECR constraints have relied so far on the hypothesis 
that protons constituted the majority of UHECRs above 10^® eV so that the step fall 
off detected at high energies could be univocally identified with the GZK one. 

Recent AUGER |184[ 1176] and Yakutsk |185] observations, however, showed hints 
of an increase of the average mass composition with rising energies np to E ^ 10^^'^ eV, 
although still with large uncertainties mainly due to the proton-air cross-section at ultra 
high energies. Hence, experimental data suggests that heavy nuclei can possibly account 
for a substantial fraction of UHECR arriving on Earth. This would be in agreement 
with fact that the initial evidence for correlations between UEHCR events and active 
galactic nuclei (AGN) — mainly Blasars — observed by AUGER |170] has not improved 
(actually has decreased) with increasing statistics |176] . 

Indeed, heavy ions are much more deviated by the extra and inter galactic magnetic 
fields due to their larger charge with respect to protons (albeit this effect is partially 
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compensated by their shorter mean free path at very high energies). The situation is 
made even more confused by the fact that these AUGER observations seem at odd with 
those of the Telescope Array (TA) based in Utah (hence in the northern hemisphere, 
while AUGER is based in Argentina) which indicate a UHECR composition dominated 
by protons [186] and some hint of correlation (the TA will need to reach about 100 
events to detect with statistic significance, at the moment about 40 vents are available) 
with large scale structures in the universe for UHECR with E > 57 EeV (albeit also 
TA does not see any significant correlation with AGN) |187] . 

Given the above described uncertainties, it is worth considering how our constraints 
on dimension 5 and 6 CPT even operators would change if the UHECR composition 
would be in the end be dominated by heavy ions. The first bad news would be the loss 
of the constraints at n = 4 on the QED sector as they are based on the proton GZK 
reaction. Actually, to be fair, this is perhaps a too harsh statement. In fact, it was 
shown in |188] that at some high energy gamma ray flux is still expected also in the case 
of mixed composition (albeit reduced with respect to the purely proton one). Hence, 
the previously discussed line of reasoning based on the absence of upper threshold for 
UHE gamma rays might still work. 

For what regards the hadronic sector, UHE nuclei suffer mainly from photo- 
disintegration losses as they propagate in the intergalactic medium. Of course, also 
photo-disintegration is a threshold process and can be strongly affected by LIV. 
According to [189] , the mean free paths of UHE nuclei are modified by LIV in such 
a way that the final UHECR spectra can show distinctive LIV features. However, a 
quantitative evaluation of the propagated spectra has not been performed yet. An 
interesting constraint was nonetheless derived in |189] assuming, for simplicity, a single 
Lorentz invariance violation parameter fj^^^ controlling Lorentz invariance violation for 
all nuclear species. A lower limit on f/^^^ was set by requiring that nuclei, which are 
stable in the Lorentz invariant case, should not undergo spontaneous decay (by emitting 
single nucleons) below energies of the order of 10^^'^ to 10^*^ eV, given that at these 
energies a substantial fraction of heavy ultra high energy cosmic ray nuclei has been 
observed. Conversely, an upper limit was given by requiring the absence of vacuum 
Cherenkov emission at the same energies. The final constraints in this case were found 
to be -3 X 10-2 < r?^^) < 4 for ^^Fe, -2 x 10-^ < n^^) < 3 x 10^^ for ^^^O and 
-7 X 10"^ < r/(^) < 1 X 10-^ for ^He, respectively. 

A final comment is deserved for the UHECR atmospheric showers reconstruction 
techniques. While preliminary, incomplete studies, hinted toward some robustness 
agains LIV modifications of the overall shower observables, there is a growing activity 
aiming at a more precise reanalysis of the problem |190[ 1191] . This is crucial to 
understand given that, for example, the absence of observation of UHE photons could 
be mimicked by anomalous shower developments which would then fool us about the 
nature of the impinging UHE particles. New results are expected soon |192] . 
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7.6. SME: Neutrino sector 

Neutrino physics can be affected in several ways by Lorentz brealcing operators. The 
most common channels for casting constraints are observations concerning the speed 
of neutrinos with respect to light, modified oscillations between neutrinos flavors and 
threshold reactions. 

Constraints from neutrino time of flights: In assessing the limit speed of neutrinos with 
respect to that of light, we have to date only a single event to rely on, the supernova 
SN1987a. This was a peculiar event which allowed to detect the almost simultaneous 
(within a few hours) arrival of electronic antineutrinos and photons. Although only few 
electronic antineutrinos at MeV energies was detected by the experiments Kamiokall, 
1MB and Baksan, it was enough to establish a constraint (Ac/c)'^^^ < 10~* |193] 
or {Ac/cY^^^ ^ 2 X 10~^ |194j by looking at the difference in arrival time between 
antineutrinos and optical photons over a baseline distance of 1.5 x 10^ ly. Further 
analyses of the time structure of the neutrino signal strengthened this constraint down to 
~ 10"^" |1951I196] . The scarcity of the detected neutrino did not allow the reconstruction 
of the full energy spectrum and of its time evolution in this sense one should probably 
consider constraints purely based on the difference in the arrival time with respect 
to photons more conservative and robust. Unfortunately adopting Ac/c < 10~*, the 
SN constraint implies very weak constraints, ^l < 10^^ and ^i, < 10^^. Note that 
observational constraints on Ac/c translate in constraints on the LIV parameter via the 
formula |105] 
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Constraints from neutrino oscillations: At odd with the previous case, we do have a 
wealth of information only about neutrino oscillations which however constraints only 
the differences among LIV coefficients of different flavors. The best constraint to date 
comes from survival of atmospheric muon neutrinos observed by the former IceCube 
detector AMANDA-II in the energy range 100 GeV to 10 TeV |197j . which searched for 
a generic LIV in the neutrino sector |198] and achieved (Ac/c)jj < 2.8 x 10~^^ at 90% 
confidence level assuming maximal mixing for some of the combinations i,j. The same 
constraint applies to the corresponding antiparticles as IceCube does not distinguish 
neutrinos from antineutrinos. IceCube is also expected to improve this constraint to 
{Ac/c)ij < 9 X 10"^^ in the next few years [199j . The lack of sidereal variations in the 
atmospheric neutrino flux also yields comparable constraints on some combinations of 
SME parameters |200] . Putting all together, it seems that no-fiavour dependent SME 
can be tolerated and many studies just assume flavour independent LIV as they starting 
base of their analysis. 
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Constraints from neutrino threshold reactions: Several threshold processes have been 
considered in the literature, most prominently the neutrino Cerekov emission u —^ '-/ u, 
the neutrino splitting u ^ uuV and the neutrino electron/positron pair production 
u — >■ z/e~e^. They are all very similar, so we shall discuss here only the latter for 
illustrational purposes and assume no LIV modification in the electron/positron sector 
as we have already seen that LIV in this sector is strongly constrained. The threshold 
energy is for arbitrary n is then 

<,„,=i^ with *,,„=e'('^y'". (80) 



The rate of this reaction was firstly computed in [104] for n = 2 but can be 
easily generated to arbitrary n [1U5] (see also |2U1] for a more general treatment and 
detailed considerations). The generic energy loss time-scale then reads (dropping purely 
numerical factors) 
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where g is the weak coupling and 9^, is Weinberg's angle. 

The observation of upward-going atmospheric neutrinos up to 400 TeV by the 
experiment IceCube implies that the free path of these particles is at least longer than 
the Earth radius implies a constraint ^ ^ 40 (taking a conservative baseline of about 
6000 Km). No effective constraint can be obtained for n = 4 LIV, however in this case 
neutrino splitting (which has the further advantage to be purely dependent on neutrinos 
LIV) could be used on the "cosmogenic" neutrino fiux which one expects from the decay 
of charged pions produced by the aforementioned GZK reaction. 

It is easy to see that the neutrino splitting should modify the spectrum of the ultra 
high energy neutrinos by suppressing the fiux at the highest energies and enhancing 
it at the lowest ones. In |202] it was shown that future experiments like ARIANNA 
|203] will achieve the required sensitivity to cast a constraint of order ^l ^ 10^^. Note 
however, that the rate for neutrino splitting computed in |202] was recently recognised 
to be underestimated by a factor 0{E/My [204] . An improved analysis can be found 
in |105] leading to an expected constraint from AUGER $^1 < 10~^ (note that this 
constraint is cast using the expected fiux which won't be able to distinguish different 
flavours) . 

However, we have just seen that experimental observations of the depth of the 
shower maximum of UHECR interactions in the atmosphere hinted at the possible 
presence of nuclei heavier than protons in UHECRs |184] . In this case, pion production 
would be suppressed at UHE and hence the UHE neutrino flux could be much smaller 
than the expectation from pure proton composition. 

A flnal comment is devoted to the so called pion decay channel tt"'" — )■ t'^yU"'' that was 
extensively explored in dealing with the recent OPERA claim |205l 1206] . In this case it 
was found that for superluminal neutrinos and unmodifled pion and muon, the process 
could be actually forbidden hence suppressing muonic neutrino pair production in the 
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CERN bean. In |2U5] it was noticed that the detection of up to 40 GeV neutrinos at 
OPERA would imply a bound ^i^^ < 10"^ Also in [2051 1206] a constraint ^l^^ < 10~^3 
was derived from atmospheric neutrinos using neutrino pair production and neutrino 
Cerenkov emission. 

While OPERA detection of superluminal neutrinos was soon recognised as flawed, 
it is worth reconsider it here as a case study. This will also illustrate some of the pitfalls 
that may be encountered when applying the above mentioned techniques to specific 
experimental situations. 

8. The Opera affaire: what was said and what could have been said 

In fall 2011 the OPERA collaboration's reported |207j and then retracted (see revised 
version of |207j ) a measurement of the time of arrival of non-oscillated z/^ over the 
path length ~ 730km from CERN to Gran Sasso with average energy of ~ 17 GeV 
which was in open conflict with Special Relativity. The u^ seemed to arrive earlier 
than light would, by an amount At = 60.7 ± 6.9 (stat.) ± 7.4 (sys.) corresponding to an 
apparent propagation velocity of Ac/c = (2.48 ± 0.28 (stat.) ± 0.30 (sys.)) x 10~^ where 
c is the low energy speed of light in vacuum, and Ac = v^ — c. After discovering a 
flaw in the initial measurement and the realisation of a new one, the OPERA result 
(Ac/c = (2.7 ± 3.1 (stat.) ± 3.4 (sys.)) x lO"*'), together with the additional results of 
Icarus |108] . and the MINOS Collaboration analysis |208] now can be seen as constraints 
on neutrino velocities in the GeV range rather than a signal of beyond standard model 
physics. 

The excitement created by the OPERA's initial report has obviously subsided for 
the greater physics community. For the quantum gravity community that focuses on 
possible experimental signatures of quantum gravity, however, technical issues were 
raised in how to analyse these types of accelerator based experiments properly. In 
particular, the detailed physics of anomalous reactions, and how they reduce the 
intensity of a particle beam from source to detector became central to the discussion. 

The most convincing theoretical objection to the initial OPERA result was 
produced by Cohen and Glashow |104] shortly after the OPERA report and involved 
just such an anomalous reaction. Cohen and Glashow used the fact that superluminal 
neutrinos should emit electron-positron pairs (see section 1716]) to argue that the OPERA 
results were not even self- consistent: any neutrino with the speed reported by OPERA 
should have lost most of its energy to pair production while it propagated from CERN 
to the detector at Gran Sasso. The maximum energy in the beam would therefore have 
dropped to be below some termination energy Et, and Cohen and Glashow showed 
that Et for the OPERA beam was less than the average ~ 17 GeV energy reported by 
OPERA. 

More specifically by integrating the energy loss rate from pair production (as 
deducible from Eq. ( IHT]) ) over a distance L and by assuming that the typical energy 
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loss length be much smaller than L one obtains 

where E is the energy on a neutrino starting with energy Eq after propagation over the 
distance L and ii^ref is the energy at which we normalize the parameter ^i,. The factor 
k = 25/448 was computed in [104j for the case n = 2, while for the general case it can be 
found in |201] . The "termination" energy Et corresponds to the energy that a neutrino 
would approach after propagation over a distance L if it started with Eq ^ Et- We 
remark here that the termination energy Et is a mildly varying function of n and of 
the energy scale -Eref- For a LIV extension of the standard model at order n = 2 (see 
the Coleman-Glashow model discussed briefly above) it was then shown that the Opera 
claim of ^^ ~ 5 X 10~^ for E^^t ~ 10 to 30 GeV implied Et ^ 12.5 GeV: a value of 
Et obviously incompatible with the observation of neutrinos up and above 40 GeV in 
Opera. Similar considerations lead to an expectation of Et — 15.0 GeV for n = 3. 

Note that this estimate was done without taking into account the similar mechanism 
of neutrino splitting. While this reaction is strictly forbidden in the n = 2 case, it is not 
so at higher orders and needs to be taken into account given that the energy loss rate of 
this process is comparable to the one for pair production loss (see e.g. |105l 12041 I201J ). 
and hence is not negligible for n > 2. In particular, this extra reaction can lead to energy 
losses without generating large numbers or electron-positron pairs as those searched for 
in an early ICARUS analysis of the Opera beam |107j . 

The physics of the Cohen-Glashow argument was correct, however the authors did 
not worry about adjusting for the finite size of the baseline. A finite baseline can be of 
the same order as the energy loss length of neutrinos undergoing pair production. This 
allows for some neutrinos to undergo only one or a few Cherenkov emissions within their 
time of flight. Therefore the most energetic neutrinos of the injection beam can still 
reach the end of the baseline with an energy larger than Et |105j . It is then necessary, 
in order to cast a robust constraint on LIV by using long baseline experiments, to run 
a full MonteCarlo simulation of the propagation of neutrinos aimed at computing the 
neutrino spectrum on arrival in the presence of this energy loss process. While this 
was not an issue for the Cohen and Glashow result, as it was one piece of a number of 
experimental and theoretical concerns about OPERA [108[ 1107] 11931 [TM] , if one wishes 
to use time of flight experiments alone to set robust constraints on neutrino LIV, the 
issue must be addressed. 

In |105] a complete analysis for the case of OPERA in the cases n = 2 and 
n = 3, taking into account neutrino pair creation and neutrino splitting (for n = 3) 
has been performed. It was there found that the propagated spectrum does indeed 
show a pronounced bump at the expected Et, but is also characterised by a high energy 
tail that extends well above Et and has an amplitude about 10% of the amplitude 
of the bump. Hence, the simple calculation of Et is not per se conclusive for casting 
constraints, although the reconstruction of the propagated spectra in [105] demonstrated 
that in the special case of OPERA the detection of neutrinos with ii^ > 40 GeV would 
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have still pointed out an incompatibility between the adopted LIV framework and the 
experimental observation. 

9. From LIV EFT to Quantum gravity phenomenology 

So far we have reviewed several constraints casted on a very conservative dynamical 
framework incorporating systematically departures from Lorentz invariance. Of course 
this has an informative value per se, given that many QG models are expected to admit 
an EFT limit well below the Planck scale. 

Nonetheless, it is also clear that the outreach of our constraints is somewhat limited 
by the very same generality of our framework and that whenever possible it should be 
now time to attempt the bolder step to apply the very same techniques developed in 
EFT to specific QG frameworks or observational/experimental evidence. 

In this sense we shall now consider a specific application of the above mentioned 
ideas with the aim to show how the past expertise could be sometime applied effectively 
to QG candidate models and observations. We shall start by first introducing an 
analogue of the mSME for the gravitational sector, i.e. the so called Einstein-Aether 
theory of gravity with its current constraints. Then we shall introduce a candidate QG 
theory which might admit a restricted version of Einstein-Aether as an IR limit and 
discuss how this can be constrained using the previous techniques, stressing in this way 
the far-fetching implications that the so-obtained constraints can have on the theory. 

9.1. Einstein-Aether theory 

The SME is constructed by coupling matter terms to non-zero LIV tensors in vacuum. 
However, when dealing with extensions of this framework to gravitational phenomena 
it is clear that leaving this tensors non-dynamical, would break general covariance. 
In this respect one can proceed along two alternative approaches, the first being 
the introduction of a suitable dynamics for the LIV tensor, the other being the the 
acceptance of an explicit breaking of the four dimensional diffeomorphism invariance of 
GR. In what follows we shall discuss the most representative models in this respect and 
they relation. 

If one choses to preserve general covariance by promoting the LIV tensors to 
dynamical fields and if restricts his/her attention to rotational invariance, then it is 
natural to generate LIV couplings by including in the action either a scalar or a timelike 
vector field that takes a vacuum expectation value. In the case of a scalar, one can use 
a shift symmetry (0(x) — )■ (pix) + (po) to construct actions for which the derivative of 
the scalar takes a non-zero value |209j . In the vector case, one simply puts a potential 
for the vector field such that the vector acquires a vev. 

We concentrate on the vector case here as it is the simplest model that allows 
for rotationally invariant Lorentz violation |210j and it is also the natural extension of 
the framework considered so far for the rotationally invariant LIV matter sector. It is 
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the most general theory for a unit timehke vector field coupled to gravity (but not to 
matter), which is second order in derivatives. 

Let us the express again the aether vector field by m" and, in analogy with what we 
did for the SME, write the most general theory for a unit timelike vector field coupled 
to gravity. If we limit ourself to second order terms in derivatives (or equivalently low 
energies) this take the form 

S = Seh + Su = —^ I d^x^/^ (R + Cu). (83) 

167rGa, J 

where £„ is given by 

Cu = -Z"%(V„M^)(V;,M^) + V{u). (84) 

The tensor Z'^g is defined as [75] 

Z'^^s = cig^'^g^s + C2S%S^s + C3S%S^^ - c^u^u^g^s , (85) 

where q, z = 1, . . . ,4 are simple coefficients of the various kinetic terms and V{u) is a 
potential term such that it generates a non-zero vev for u°'. Note the indicial symmetry 
Z r = Z g . An additional term, RabU°'u^ is a combination of the above terms when 
integrated by parts, and hence is not explicitly included here. 

The potential in f lM|) is normally fixed by the requirement to remove the ghost 
excitation associated to the one of the vector components which will necessarily acquire 
a wrong sign in the kinetic term. The choice V{u) = A(m^ + 1), where A is a Lagrange 
multiplier, fixes the norm of m" and removes the ghost excitation (c.f. the discussions 
in [71] and [211] ). It is this choice of the potential that it is normally associated to the 
so called "Einstein- Aether theory" [210] which we shall consider here. 

9.1.1. Constraints on Einstein- Aether gravity: For what regards the constraints on 
Einstein-Aether gravity, they can be divided in those on the aether kinetic terms and 
those on the aether-matter couplings. Given the the latters can be reduced to mSME 
constraints, we shall here discuss specifically only the first kind (one can find a discussion 
of the second kind in e.g. [120] ). 

Before of doing so, we can note however that the couplings between aether and 
the standard model field content while being the same as those discussed before in this 
review, they do have new features as the aether field has now dynamics, so there can be 
position dependent violations of Lorentz symmetry. Interestingly, some of the couplings 
to matter that are unobservable for a single fermion field can have relevant effects when 
the aether varies. For example, the —aunip'-f'^ip term in the mSME could be removed 
by making a phase change for the fermion. However, once m° is dynamical and varies 
with position, only a single component of the term can actually be removed by a phase 
change [212] . This leads to a new type of term which requires gravitational/position 
dependent tests in the matter sector [212] . 

In order to constrain instead aether kinetic terms one can adopt two different 
approaches. One consists of course in adopting, as normally done any modified gravity 
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theory, a PPN analysis that allows to compare the theory with observations e.g. solar 
system constraints. Alternatively, one can use the fact that the theory predicts extra 
degrees of freedom in the gravity sector (there are naively four, but the unit constraint 
removes one) and that these excitations strongly couple to the metric via the unit 
constraint (for a more detailed discussion see [751 1120] ). 

Constraints from PPN analysis: All the PPN parameters vanish except for ai, 0:2 which 
describe preferred frame effects, ai and 02 were calculated in |213] 

-8(4 + C1C4) 
"^ = 2c, - c? + ci ^^'^ 

«! (Ci + 2C3 - C4) (2Ci + 3C2 + C3 + C4) 

02 = -, 77 r ■ [o7) 

2 (Ci +C2 + C3)(2-Ci -C4) 

Current constraints are ai < 10~^ and 02 < 4 x 10^^ |214] and so from a PPN analysis 
alone there is still a large 2-d region of parameter space that remains consistent with 
available tests of GR. 

Constraints from gravity- aether wave modes: The combined aether-metric modes 
consist of the two usual transverse traceless graviton modes, a vector mode, and a scalar 
mode |215j . The speeds of each of the modes can differ from the speed of light. Hence if 
the speeds are less than unity, high energy cosmic rays will emit vacuum gravitational 
Cerenkov radiation |216] . If we denote the speeds of the spin-2, spin-1 and spin-0 modes 
by S2, Si, So then we have |215j 

sl = {l-cr~ cs)-' (88) 

„2 _ 2ci -cl + cl ^gg^ 



2(ci + C4)(1 -Ci -C3) 

(Ci + C2 + C3)(2-Ci -C4 



"'^ (Ci + C4)(l-Ci-C3)(2 + 3C2 + Ci+C3)' ^^^^ 

The requirement that all these speeds are greater than unity therefore puts constraints 
on a combination of the Cj coefficients. However, even after imposing all of the above 
constraints there is still a large region of parameter space allowed. Indeed, the PPN 
and gravitational Cerenkov constraints are all satisfied provided quite lose conditions 
on the model coefficients are satisfied [751 I120J . Hence, the gravitational sector is only 
minimally constrained compared to aether-matter couplings. 

9.2. Hofava-Lifshitz gravity 

The underlying idea of the Hofava-Lifshitz (HL) gravity (see e.g. [217] for a review) 
is to achieve power-counting renormalizability by modifying the graviton propagator in 
the ultraviolet by adding to the action terms containing higher order spatial derivatives 
of the metric, but not higher order time derivatives, so to preserve unitarity. This 
procedure naturally leads to a space-time foliation into spacelike surfaces, labeled by 
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the t coordinate and with x' being the coordinates on each surface. The resulting theory 
is then invariant only under the reduced set of diffeomorphisms that leave this foliation 
intact, t — )■ t{t) and x* — )■ x*(t, x*). 

It was shown that power counting renormalizability requires the action to includes 
terms with at least 6 spatial derivatives in 4 dimensions [761 1218] . Of course, all 
lower order operators compatible with the symmetry of the theory are expected to 
be generated by radiative corrections, so the most general action takes the form 



Shl = ^I dtd'x NVh (^L, + ^L, + ^Le) 



(91) 



where h is the determinant of the induced metric hij on the spacelike hypersurfaces, 

L2 = KijK'^ - XK^ + ^^''^R + r]aia' , (92) 

where K is the trace of the extrinsic curvature Kij, ^^^R is the Ricci scalar of hij, N is 
the lapse function, and Oj = di In A^. L4 and Lq denote a collection of 4th and 6th order 
operators respectively and M^ is the scale that suppresses these operators which does 
not coincide a priori with Mpi. 

It is perhaps tempting to call M^ the Lorentz breaking scale, but the theory exhibits 
Lorentz violations (LIV) at all scales, as L2 already contains LIV operators. These 
Infrared (IR) Lorentz violations are controlled by three dimensionless parameters that 
take the values A = 1, .^ = 1 and r/ = in General Relativity (GR). While, C, can be 
set to 1 by a suitable coordinate rescaling, it is presently unclear if the running of the 
remaining two parameters will converge on the GR values in the IR. Nonetheless, it 
seems that the the theory could still be viable and consistent for suitable choices of the 
dimensionless parameters A, rj which admits the GR values as extremal limits of the 
allowed range j219] . 

Action f l9T|l does present, however, the unappealing feature to contain (in L4 and 
Lq) a very large number, O(IO^), of operators and independent coupling parameters. 
In remedy of this situation, restrictions to the theory have been proposed which would 
limit the proliferation of independent couplings. We shall not deal with such restrictions 
here (but see e.g [2171 1220] for a concise review) as they will not be determinant for the 
phenomenological discussion on HL that we shall present later on in this review. 

9.3. Relation between Einstein-Aether and Hofava-Lifshitz gravity 

One interesting, and at the same time problematic, feature of HL gravity is the presence 
of a new propagating scalar mode associated with the reduced diffeomorphism invariance 
of the theory with respect to GR. However, if one chooses to restore diffeomorphism 
invariance, then this mode manifests as a foliation-defining scalar field [79]. This field 
allows also to make manifest the relation between the L2 Lagrangian and the Einstein- 
Aether one. 

In fact, it was shown in [79] (but see also [221j ) that the Einstein-Aether theory 
is equivalent to the the infrared limit of HL gravity if the aether is assumed to be 
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hypersurface orthogonal before the variation. More precisely, hypersurface orthogonality 
can be imposed through the local condition 

where T is a scalar field that defines a foliation. Choosing T as the time coordinate one 
selects the preferred foliation of HL gravity, and the action (l83l) reduces to the action 
of the infrared limit of HL gravity, whose Lagrangian we denoted as L2 in eq. flOTj) . 
The details of the equivalence of the equations of motions and the correspondence of 
the parameters of the two theories is discussed in detail in |222l I223J . 

As a concluding remark let us stress that the Einstein- Aether theory can be seen as 
a sort of gravity sector version of the rotational invariant mSME. However, one should 
expect to be able to supplement also Einstein-Aether gravity with suitable higher order 
operators. This has not been done yet, but the existence of HL gravity and its relation 
in the IR with Einstein- Aether theory seems to suggest that such an extension should be 
viable. For example, a possible extension of HL gravity to a theory invariant under the 
full group of four-dimensional diffeomorphisms have been proposed |224j . This kind of 
extensions could then be seen as possible generalisations of the Einstein- Aether theory 
beyond the IR limit. 

9.4- Constraints on Hofava gravity 

Coming to constraints on HL gravity theory, it should be obvious that given the relation 
in the IR between hypersuface orthogonal Einstein- Aether and Hof ava-Lifshiftz gravity, 
the previously presented constraints on Einstein-Aether gravity can in principle be 
related to constraints for the latter theory. 

Looking then at the UV complete the theory from the point of view of QG 
phenomenology, it then interesting to know the available constraints on the Lorentz 
breaking scale of the theory M^. Remarkably, this scale happens to be bounded both 
below and above. Indeed, for the theory to preserve power counting renormalizablity 
and be at the same time compatible with current observations (microgravity experiments 
and solar system tests) one has to require 0(1) meV < M^ < 10^^ GeV [2T9l [2T71IT57] . 
a quite broad opportunity window for the theory. 

We have seen however, that radiative corrections will always allow for LIV operators 
to percolate from a SM sector to the other, gravity being no exception. In this sense it 
can be generically expected for the above theory to induce Lorentz breaking operators 
in the matter sector at all orders. Let us assume here that no LIV is present in the 
matter sector at tree level and that again some protective mechanism will prevent the 
percolation of the Lorentz breaking terms to the lowest order (mass dimension 3 and 4) 
operators of the matter sector (we shall come back later on this point). Also, one can 
assume that the CPT and Parity (P) invariance of the gravitational action is preserved in 
the matter sector. Indeed if no CPT and P odd operators are present in the matter sector 
at the tree level one would not expect them to be generated via radiative corrections 
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induced by the gravitational, CPT and P even, terms. This assumption forbids hehcity 
dependent terms and allows only even power of the momentum in the matter dispersion 
relation. Hence matter and anti-matter are expected to share the same dispersion 
relation which within this framework one can then expect to be 

From a logic point of view, there are now two options: (a) M* = Mliv, i.e. M* is a 
universal scale; (b) M^ <^ Mliv (a^s we have seen in this review current phenomenological 
constraints already rule out the case Muv <^ M^). Clearly, option (b) requires some 
mechanism which suppresses the percolation of LIVs in the matter sector even to higher 
order operators. We already discussed such a mechanism, the so called "gravitational 
confinement" [131] . Here, reporting the work done in [157j . we shall conclude the 
necessity for such mechanisms. 

For doing so the strategy adopted in |157j consisted in focussing on option (a) and 
demonstrating that, in this case, matter LIV constraints imply M* = Mliv > 10^^ GeV, 
thus closing the available window for M^. Of course this could be easily achieved using 
the constraints on QED for modified dispersion relation of order n = A using UHECR. 
However, we have seen that these constrains are somewhat questionable nowadays while 
further evidence about the nature of the highest energy particles is awaited. 

Hence, in |157j . it was suggested to use the more robust observation of the 
synchrotron radiation from the Crab Nebula. A first estimate of the strength of the 
available constraint can be derived along the lines described before. In particular for a 
negative t] one can run the argument that the maximal obtainable frequency ( 1761) should 
not be smaller than the observed maximal one. 

The maximal observed frequency in the CN synchrotron spectrum, Uohs ~ 0.1 GeV, 
for n = A one can then analytically derive a constraint r] > —10^ (assuming B ~ 300 /iG 
and Mliv = Mpi), which would correspond to Mliv > 3 x 10^® GeV. 

While promising, this is not a double sided constraint. In order to be sensitive to 
the full range of the rj parameters and take into account competing LIV and LI effects 
a much deeper analysis is needed. As we saw, this was performed in |156j where the 
possible LIV induced modifications to the standard Fermi mechanism (which is thought 
to be responsible for the formation of the spectrum of energetic electrons in the CN) 
were considered and the synchrotron spectrum of the CN was recomputed taking into 
account all the new, LIV induced, phenomena for n = 3 LIV. This was done in [157] for 
n = 4 and the specific modified dispersion relation Eq. (IMj) . 

The free parameters of the model (electron/positron density and spectrum and 
magnetic field strength) were fixed in order to reproduce the low energy part of the 
spectrum, which is not affected by LIV This allows to reconstruct how LIV affects the 
higher energy part of the spectrum {E >100 keV) and to use a x^ statistics to measure 
when deviations from the observed spectrum due to LIV become unacceptably large. 
By considering the offset from the minimum of the reduced x^ exclusion limits at 90%, 
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95% and 99% Confidence Level (CL), according to |225] it was shown in |157] that mass 
scales Mliv < 2 x 10^^ GeV are excluded at 95% CL. 

So one can conclude that the application of the previously described methods to 
concrete QG model leads to a couple of important lesson. First of all, one should 
not think a priori that some methods are effective only for some order of LIV EFT. 
The synchrotron constraint we considered here is apparently ineffective for n = 4 and 
Mliv = ^Ph however it is not obvious that any QG model should send the LIV scale 
equal to the Planck one. Secondly, the constraint obtained has broader implication for 
the model that could not be foreseen in a purely EFT approach. In particular it implies 
necessarily M^, <^ Mliv and hence that a protective mechanism beyond the lowest order 
operator should be foreseen in the matter sector of Hofava gravity. As we discussed 
previously it has been shown that dimension 4 operators of matter can be efficiently 
screened from LIV this way in HL models |131] and this is expected to be the case for 
higher order operators as well. The investigation we just reported can then be taken 
as a strong indication that such a mechanism should be considered the main avenue for 
the making this particular QG proposal viable. 

10. Beyond EFT: other frameworks 

Specifying which dynamical framework is employed is crucial when discussing the 
phenomenology of Lorentz violations. In this review we have focussed on the the 
most conservative framework, effective field theory. However, this is not the only 
one conceivable and, in fact, there are reasonable arguments from holography that 
a quantum gravity theory should not necessarily be a local field theory in the UV 
(c.f. the discussion in |226] ). Hence deviation from standard Lorentz invariance may 
enter into low energy physics in novel ways. Given that the EFT approach is nothing 
more than a highly reasonable, but rather arbitrary "assumption" , it is worth studying 
and constraining additional models, given that they may evade the majority of the 
constraints discussed in this review. 

10.1. D-brane models 

A class of models showing modified dispersion relations of the general form Eq. (l52il 
was derived from the Liouville string approach to quantum space-time |227[ [50l 152] . 
These models motivate corrections to the usual relativistic dispersion relations that are 
first order in the particle energies and that correspond to a vacuum refractive index 
?7 = 1 — (E'/Mpi)", where a = 1. Models with quadratic dependences of the vacuum 
refractive index on energy: a = 2 have also been considered [67] . 

Most importantly, the D-particle realisation of the Liouville string approach predicts 
that only gauge bosons such as photons, and not charged matter particles such as 
electrons, might have QG-modified dispersion relations. This occurs since excitations 
which are charged under the gauge group are represented by open strings with their ends 
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attached to the D-brane |228] , and only neutral excitations are allowed to propagate in 
the bulk space transverse to the brane |229] . Thus, if we consider photons and electrons, 
in this model the parameter r) is forced to be null, whereas ^ is free to vary. Noticeably, 
the theory is also by construction CPT even, implying that vacuum is not birefringent 
for photons {C,+ = ^-)- As such these models are very difficult to constraint. However, 
one of the possible incarnations of this model with a = 1 was constrained in |230] by 
using again the absence of upper threshold for secondary photons of the GZK reaction. 
The obtained constraint on the photon LIV parameter is of order 0(10^^^) however 
alternative spacetime foam models for which this constraint does not apply can be 
envisaged [231] . 

10.2. New Relativity Theories 

An alternative route with respect to explicit Lorentz breaking pursued in the extant 
literature consists in trying to define a new relativity group which would allow for an 
extra invariant scale (generally the Planck length) beyond the speed of light. In this 
direction is then worth exploring generalisation of the von Ignatowski derivation which 
would keep the relativity principle while relaxing one or more of the above discussed 
postulates. 

For example, relaxing the space isotropy postulate leads to the so-called Very 
Special Relativity framework [232], which was later on understood to be associated 
to a Finsler-type geometry [2331 12341 I235J . In this example, however, the generators 
of the new relativity group number fewer than the usual ten associated with Poincare 
invariance. Specifically, there is an explicit breaking of the 0(3) group associated with 
rotational invariance. 

Much research has instead focus on constructing alternative relativity groups 
with the same number of generators as special relativity. Currently, we know of no 
such generalisation in coordinate space within the standard commutative geometry 
framework. However, it has been suggested that, at least in momentum space, such a 
generalization is possible, and it was termed "doubly" or "deformed" (to stress the fact 
that it still has 10 generators) special relativity (DSR) [871 [SHI [89l [90] . Such momentum 
space construction have been linked to possible Lie-type non commutative geometry 
(mainly the so called kappa-Minkowski one |236[ 1237] ) or to curved momentum spaces 
in higher dimensions [238] . Even though DSR aims at consistently including dynamics, 
a complete formulation capable of doing so is still missing, and present attempts face 
major problems. 

Also to overcome some of these conceptual problems was recently advanced what 
we might call a "spin-off" of DSR named Relative Locality. This framework is based on 
the idea that the invariant arena for classical physics is a curved momentum space rather 
than spacetime (the latter being a derived concept) [HH [921 [931 IHH [95]. Interestingly, 
while this seems to allow an extension of the relativity principle with an invariant energy 
scale, it does so at the price to renounce to an observer-independent concept of locality. 
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There is some evidence that this might be generic feature of any alternative relativity 
group of this sort [93], suggesting locality violations in EFT as a possible new avenue 
of exploration in QG phenomenology (see also discussion on non-locality in QG below). 
In summary, DSR and Relative Locality are still a subject of active research and 
debate (see e.g. |239l 12401 12411 1242] ): nonetheless, they are reaching just now the level 
of maturity required for casting constraints (see e.g. |243] ). 

11. Discussion and perspectives 

We summarize the current status of the constraints for the LIV SME (rotational 
invariant) in Table El Of course at first sight this might seem a quite satisfactory 



Order 


photon 


e-/e+ 


Protrons 


Neutrinos'* 


n=2 
n=3 

n=4 


N.A. 

0(10-16) (GRB) 

0(10-8) (CR) 


0(10-16) 
0(10-16) (CR) 
0(10-8) (CR) 


0(10-20) (CR) 
0(10-14) (CR) 
0(10-6) (CR) 


0(10-8^10-10) 
0(40) 
0(10-^)* (CR) 



Table 2. Summary of typical strengths of the available constrains on the SME at 
different n orders for rotational invariant, neutrino flavour independent LIV operators. 
GRB^ganima rays burst, CR^cosmic rays. " From neutrino oscillations we have 
constraints on the difference of LIV coefficients of different flavors up to 0(10-^8) on 
dim 4, 0(10-8) and expected up to 0(10-1"^) on dim 5 (ICE3), expected up to 0(10-'') 
on dim 6 op. * Expected constraint from future experiments. 



state of the art, so much so that one might ask if we haven't tests Lorentz violations 
enough and should now move one towards new phenomenology. As usual, the answer is 
not a sharp one. Let us further elaborate on this point. 



11.1. Uncertainties on n = A constraints 

Let's first stick to tests of violation of Lorentz invariance in the SME. Here, as we 



discussed at length in section 17.51 the main open issue is provided by the lasting 
uncertainty about the UHECR composition and heck the actual observation of the GZK 
cutoff. In this respect the following comment is in order. The observational picture is yes 
confused but not hopeless. The issue will be probably settled in a few years and some 
experiments like TA seem still to provide a more conservative picture than AUGER. 
As a matter of fact, it still seems more probable that in a few years, the constraints 
we presented in this review based on the GZK cutoff will be confirmed or strengthened 
rather than disproved. 

It would be however unfair to play down the present uncertainties in UHECR 
physics and place this constraints at the same level of robustness e.g. of those cast 
at order n = 3 by using the synchrotron radiation from the Crab nebula. In this 
sense any UHECR-independent constraint (by which we mean limiting the relevant LIV 
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parameters to be less than one, being this their natural value within the SME) on the 
dimension 6 operators of the SME would be more than welcomed. As we have seen in 
section 19. 4[ this is possible within specific QG models, but so far only when the Lorentz 
breaking scale is theoretically already set to be smaller than MpL. 

11.2. The naturalness of Lorentz violations 

Another open issue is of course the naturalness problem. Lacking so far evidence for 
new physics at intermediate scales from the Higgs till the Planck scale one might wonder 
if the only hope for a working Lorentz breaking theory is gravitational confinement (see 
section 15. 2p with its consequent issue of large LIV in the gravitational sector. It is 
too early to say, but of course investigations on specific models like the one presented 
in section 19.41 seems to be strongly suggestive towards this direction. It is probably 
too early to say, but it is clear that the naturalness problem is probably the most 
pressing theoretical challenge LIV models are facing today and almost a selection tool 
for candidate theories which admit a low energy EFT description. Furthermore, we 
stress again that any emergent gravity scenario will have to be predictive in this sense 
as relativistic behaviour of its fundamental constituents does not seem enough, per se, 
to guarantee an exact relativistic emergent system |12j . 

11.3. Towards an authentic QG phenomenology? 

Perhaps a final comment is due to the present state of the art of the field. Lorentz 
breaking phenomenology has been a remarkable success, a community effort which has 
built (in a bit over a decade) a wealth of knowledge, methods and constraints that are 
now at our disposal for efficiently testing candidate theories of QG once their low energy 
limit is known. However, this state of affair is not and cannot be the end of the story as 
EFT phenomenology is not yet a true quantum gravity phenomenology. Examples like 
the one discussed in section 19.41 should have convinced the reader that the implications 
of the phenomenological constraints that can be casted in full fledged QG models are 
more far reaching than those derived in the purely phenomenological EFT approach to 
LIV. We should further pursue this line of research and as ask more forcibly now to our 
QG models for testable low energy predictions. 

There are probably many more phenomenological consequences of QG beyond and 
apart LIV and maybe some of them will go beyond the realm of local EFT. For example, 
causal sets models seems to entail non-locality as a counterweight to the requirement 
of a Lorentz invariant discretisation of spacetime [244j . as such the local nature of the 
low energy world should be in these models emergent and not exact. Can we test such 
non-localities? in which contexts they could be probed? Let us note again, as a side 
remark, that a link between some form of non-locality and extension of special relativity 
in situations which seem to violate the homogeneity axiom of von Ignatovski derivation 
have similarly been suggested [93] and probably deserve further investigations. 
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In conclusion, there is probably a long way to go and the challenges ahead are 
not just for the QG phenomenology community but also for those working on QG 
models. The very existence of the above described constraints on Planck scale new 
physics requires a new effort from these communities to search more actively for new 
windows for testing quantum gravity. The best has yet to come. 
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